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HAUPTAUFSÄTZE 


Exakte Lösungen der Differentialgleichungen einer 
 adiabatischen Gasströmung. 
Von Friedrich Ringleb (Messerschmitt A.G.) in Augsburg. 


vn Lösungen der Differentialgleichungen einer adiabatischen Gasströmung liegen bisher 
nur vereinzelt vor‘). Überdies sind die vorhandenen Lösungen nicht ausreichend ge- 
wesen, genügend Einsicht in die Verhältnisse zu bringen, welche beim Durchgang der Strö- 
mungsgeschwindigkeit eines Gases durch die Schallgeschwindigkeit und im Überschallgebiet 
herrschen. Insbesondere dürften Begriff und Ursache des Verdichtungsstoßes nicht hinreichend 
geklärt sein. 

Im folgenden werden die Differentialgleichungen der adiabatischen Gasströmung in ein- 
facher Weise exakt integriert. Dabei gelingt es, besonders einfache Beispiele für Gasströ- 
mungen zu ermitteln. Als hauptsächlichste Ergebnisse der folgenden Untersuchungen seien 
genannt: 


1. Bei einer Strömung über eine gewölbte Fläche ist die Schallgeschwindigkeit im all- 
gemeinen kein kritischer Zustand ?). 

2. Kritisch sind nur Strömungsstöße. Unter einem Strömungsstoß verstehen 
wir dabei eine Stelle, an der die Beschleunigung der adiabatischen Strö- 
mung unendlich wird. Im allgemeinen bildet die Gesamtheit der Stoßstellen eine Kurve, 
die wir kurz Stoßlinie nennen. 

3. Längs einer Stoßlinie besitzen die Stromlinien Rückkehrpunkte (Spitzen). 


4. Die Geschwindigkeitskomponente senkrecht zur Stoßlinie ist gleich der Schall- 
geschwindigkeit. 


Dieser theoretische, gesamte Strömungsverlauf ist natürlich praktisch unmöglich. Die 
Strömung wird in Wirklichkeit im allgemeinen schon vor der Stoßlinie ihren adiabatischen 
Charakter verlassen, also z. B. mit einem Verdichtungsstoß in eine mögliche Unterschall- 
strömung übergehen, oder sie wird sich selbst zerstören. Als praktisch realisierbar sehen wir 
aber jeden strömungsstoßfreien Teil einer theoretisch ermittelten adiabatischen Strömung an. 


1) Vgl. A. Busemann: Gasdynamik, Handbuch der Experimentalphysik Bd. 4, Teil 1, Leipzig 1931, S. 443. 
— J. Ackeret: Gasdynamik, Handbuch der Physik Bd.7, Berlin 1927, S. 288. — L. Prandtl: Allgemeine Überlegungen 
über die Strömung zusammendrückbarer Flüssigkeiten. Convegno di science fisiche, matematiche e naturali. Atti 
dei convegni 5 (1936) S. 169 bis 197 und die dort zusammengestellte Literatur. — W.Tollmien:Zum Übergang von Unterschall- 
in Überschallströmungen. Z. angew. Math. Mech. Bd. 17 (1937), S. 117 bis 136. — A. Busemann: Hodographenmethode 
der Gasdynamik. Z. angew. Math. Mech. Bd. 17 (1937), S. 73 bis 79. — C. Ferrari: Campo di moto in fluidi com- 
pressibili. Aerotecnica Bd. 18 (1988), S. 400 bis 411. 

2) Darauf hat schon Tollmien hingewiesen a. a. O. 
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1. Die Differentialgleichungen einer adiabatischen Gasströmung. Wir bezeichnen mit x, y 
die rechtwinkligen Koordinaten eines Punktes in der Ebene, in welcher die reibungslose, 
stationäre adiabatische Gasströmung erfolgen soll. «, ® seien die rechtwinkligen Koordinaten 
des Geschwindigkeitsvektors der Strömung, w, ö seine Polarkoordinaten. w ist also der ab- 
solute Betrag der Geschwindigkeit, 9 die Strömungsrichtung und es gilt 


cos», v=mwsin®. 


Die Schallgeschwindigkeit des ruhenden Gases wird mit e,, die örtliche Schallgeschwindigkeit 
bei der Strömungsgeschwindigkeit w mit ce bezeichnet. Dann ist 


wobei x für zweiatomige Gase (Luft) gleich 1.4 zu setzen ist. Ferner sei o, die Dichte des 
ruhenden Gases, o die Dichte des mit der Geschwindigkeit w strömenden Gases. Es gilt 


1 
% 2 
Zwischen dem Druck p, des ruhenden, dem Druck p des strömenden Gases, seiner Geschwin- 
digkeit w und ce, besteht schließlich die Gleichung 


Die Differentialgleichungen der Strömung lauten’): - 


wobei die erste Gleichung die Wirbelfreiheit, die zweite Gleichung die Quellen- und Senken- 
freiheit (Kontinuität) der Strömung bedingt. Wird das strömende Medium inkompressibel an- 
genommen, so gehen wegen o=o, die Gl. (4) und (5) in die Cauchy-Riemannschen 
Differentialgleichungen für « und — v über, deren kompressible Verallgemeinerungen sie sind. 


Ebenso wie im inkompressiblen Fall folgt aus (4) und (5) die Existenz zweier Funktionen 
D=P(x,y), Y=Y(x,y) derart, daß 


u—®,, 


gilt. Wir nennen ® das Geschwindigkeits-, Y das Strömungspotential. 


Die Kurven /=const sind die Stromlinien. Denn für diese ist die Tangentenrichtung 


de 
Andererseits folgt für die Kurven Y=const durch Differentiation 


Nach Gl. (7) ist aber die rechte Seite gleich — 


Elimination von « und v aus den Gl. (6) und (7) liefert die Bedingungsgleichungen 
fürd® und 


d,=-7, . . . . . . . . (8), 


aus welchen sich die Gleichungen 


3) Die Indizes x und y bedeuten ebenso wie a 50 Differentiation nach x bzw. y. 
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ableiten lassen *). In diesen kann man noch «, v» und e mit Hilfe der Gl. (6), (7) und (1) durch 
die Ableitungen von ® bzw. Y ausdrücken. Für eine mathematisch exakte Lösung sind diese 
Gleichungen sehr kompliziert und wenig brauchbar. Dagegen kann man, wie L. Prandtl 
dies getan hat, die « und » und damit auch ce in den Koeffizienten der Gleichungen im Kleinen 
konstant annehmen und dann die entstehenden linearen Gleichungen mit konstanten Koeffi- 
zienten in bekannter Weise lösen. 


Die Bedingungsgleichungen für ® und Y lassen sich aber auch im 
Großen und zwar exakt in lineare überführen, indem man w und 9 als un- 


abhängige Variable an Stelle von x und y einführt°). Führt man nämlich in den 
Gleichungen _ 


dP=ude+rdy, 
vdetudy) 
00 


für u und » die Polarkoordinaten w, ® ein und löst die Gleichungen nach dx und dy auf, 
so erhält man 


o w 


Für eine vorliegende Strömung kennt man an jeder Stelle x, y die Geschwindigkeit w und 
ihre Richtung d. Es sind also w und # Funktionen von &, y und umgekehrt. Denken wir 
uns x, y als Funktionen von w und Ö, so werden auch ® und Y Funktionen von w und Ö. 
Aus den letzten Gleichungen erhält man alsdann die Beziehungen: 


cos d 0, sind 


cos d o, sin d 
(11). 
sin d 0, 
sin d 0, cos" 
Es müssen nun die Bedingungen 
Yo w 


erfüllt sein. Differenziert man also die erste und dritte Gleichung von (11) nach 9 und die 
übrigen nach w und setzt die erhaltenen Ausdrücke in den letzten Gleichungen ein, so er- 
geben sich, wenn man bei der Rechnung die Beziehung 


d 
beachtet, für ® und Y die beiden Bedingungen: 
w’\ 1 


4) W genügt nicht etwa derselben Differentialgleichung wie ®, denn u und ® enthalten ja noch die abhängige 
Variable auf verschiedene Weise. 


5) Dieser Übergang findet sich bereits bei P. Molenbroek: Über einige Bewegungen eines Gases bei Annahme 
eines Geschwindigkeitspotentials. Arch. d. Mathem. u. Phys., Grunert Hoppe 18%, Reihe 2, Bd.9, S. 157. — A. Tscha- 
pligiu: Wiss. Ann. d. Univ. Moskau, Math. Phys., Bd. 21 (1904), S.1 bis 121 (russisch). Vgl. dazu D. Riabouechinski: 
Comptes rendus Acad. Sei., Paris Bd. 194 (1932), S. 1215; B. Demtschenko: ebenda Bd. 194 (1932), S. 1218 und 1720, 
sowie Publ. Math. Univ. Belgrade Bd. 2 (1933), S. 85. — A. Steichen: Beiträge zur Theorie der zweidimensionalen 
B:wegungsvorgänge in einem Gase, das mit Überschallgeschwindigkeit strömt. Diss. Göttingen 1909. 
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Diese Gleichungen entsprechen völlig den Gl. (8) für ® und Y als Funktionen von x, y. 
Differenziert man schließlich noch die erste dieser Gleichungen nach 9, die zweite nach w. 
so folgt für Y die Differentialgleichung‘) 


w? 


Für die Integration hat man noch e durch w auszudrücken und erhält 
wobei zur Abkürzung 
”—3 
gesetzt ist. 
Zur Ermittelung von Strömungen hat man folgenden Weg einzuschlagen: 
Man sucht eine Lösung Y der linearen Differentialgleichung (15), bildet 7, und 3 und mit 
Benutzung von (13) auch ®,. und ®3. Die Gl. (11) liefern dann x,., 29, Yw, 9a als Funktionen 
von w und 9. Aus ihnen erhält man x und 4 selbst durch einfache Quadraturen. Man hat 
somit x und y als Funktionen von w und ®: 
y=y(w,d). 
Eine Stromlinie befriedigt, wie auf S. 186 unten gezeigt wurde, die Bedingung 
Y(wd)=k, 


wo k eine Konstante ist. Berechnet man hieraus z. B. 9 als Funktion von w und setzt oben 
ein, so ergibt sich für jede Stromlinie eine Parameterdarstellung 


2—=f(w, k), y=g(w,k). 


Eine Folge der Linearität der Gl. (14) bzw. (15) ist die Möglichkeit, Lö- 
sungen zu überlagern. Sind nämlich 7, und Y, Lösungen von (15) und ®,, ®, ver- 
möge (13) zugehörige Geschwindigkeitspotentiale, so ist auch €, %, +C, Y, (C,, €, Konstante) 
eine Lösung von (15) und €, ®,+(,®, das zugehörige Geschwindigkeitspotential. Sind dann 
2,9, zu ®D,,Y, und x,,9y, zu ®,, Y, gehörige Lösungen von (11), so sind 

=0,%,+0,%,, 
Funktionen von w, d, welche zusammen mit der Gleichung 


+0,%,=k 


das System der Stromlinien einer neuen adiabatischen Gasströmung liefern, wenn k alle mög- 
lichen konstanten Werte annimmt. 


2. Quelle, Wirbel und Wirbelquelle. Begriff des Strömungsstoßes. Wir wenden das ge- 
schilderte Verfahren zunächst auf einige besonders einfache Fälle an. 


Es sollen erstens diejenigen Strömungen bestimmt werden, deren sämt- 
liche Stromlinien einem Geradenbüschel angehören. 


Für jede Stromlinie ist 9 konstant. Da die Stromlinien andererseits der Bedingung 
Y= const genügen, kann also Ylediglich eine Funktion von # allein sein. Aus der Differential- 
gleichung (15) für Y folgt dann sofort 


dw 
und 
YV=(09+0' (0, C' konstant). 
Weiter ist 


Die Gl. (13) liefern daher 


—0. 


6) Vgl. A. Steichen, a. a. 0. 


2. . Math. h. . A . 
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Die Gl. (11) ergeben 


0%, m’\ cos d 0, sind 
o w o w 
w o w 


Aus der zweiten und vierten Gleichung erhält man durch Integration nach # 


0, 08 0, sin d 
= w 
o w +fw), o w 


+9(m), 


wobei die noch unbekannten Funktionen f(w) und g (mw) die Integrationskonstante vertreten. 
Sie werden mit Hilfe der ersten und dritten Gleichung bestimmt. Setzt man nämlich die ge- 
fundenen Ausdrücke für « und y in diesen Gleichungen ein, so findet man 


also 
0, 0, sin d R 
z=0=-—— +K 
x w + 17 Y o +K,, 


wo K, und K, Konstante bedeuten. Diese Gleichungen geben für die gesuchte Strömung den 
Zusammenhang zwischen den Koordinaten w, # der Geschwindigkeit und dem Ort x, y an. 
Die Stromlinien selbst genügen der Bedingung 


(k konstant). Wir erhalten somit für die Stromlinien die Parameterdarstellung 


(C” konstant) oder 


0, cos k o, sink 
K y=C- K 
oe w o w 


mit dem Parameter w. Ohne wesentliche Beschränkung der Allgemeinheit setzen wir noch 


C=c, K=K,=0 
und erhalten 


-sin k. 
o 


Die Kurven konstanter Geschwindigkeit sind die Kreise 


2 
€, 
ow) 


In Bild 1 ist das Impulsverhältnis eT & T 
co7o 40 T 


graphisch dargestellt. Für im 


| 

1,095’ 
also wenn w gleich der örtlichen Schall- 45 
geschwindigkeit wird, erreicht die 
Kurve ihr Maximum. Wenn nun w 
alle Werte von O bis zur Maximal- 
geschwindigkeit 2,236 e, durchläuft, so 
nimmt der Radius des obigen Kreises 
om 
ab und wächst dann wieder von diesem Werte bis ©. Im Gegensatz zur Quelle (bzw. Senke) 
eines inkompressiblen Mediums ist diese jetzt nicht mehr punktförmig, sondern ein Kreis 
(Grenzkreis”?)), dessen Radius gleich dem reziproken Werte des Impulsverhältnisses ist. Die 
Strömung erfährt an dem Grenzkreis einen Stoß, denn man hat, wenn ds das 
Bogenelement der Stromlinie bedeutet, 


20 6236 


von © bis zum Minimum 


\2 „2\2 
(1-5) 


e? 


dw?, 


?) Soweit ist die Gestalt der kompressiblen Quelle bzw. Senke bekannt. 


| 
w 
| | 
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also 


Andererseits ist 


wo b die Beschleunigung in der Bahnrichtung bedeutet. Für den Grenzkreis gilt also 
w=c, 


d. h. aber, der Grenzkreis ist eine Stoßlinie°). Die Strömung würde nun mit der gleichen 
Geschwindigkeit in entgegengesetzter Richtung zurückprallen und ihren Verlauf mit steigender 
Geschwindigkeit fortsetzen, bis sie im Unendlichen die Maximalgeschwindigkeit erreichen würde 
(oder umgekehrt), wenn sie dazu die Möglichkeit besäße. Das ist aber offenbar nicht der 
Fall, da die Strömung stationär vorausgesetzt ist und die reflektierte Strömung gegen die 
nachfolgende noch nicht reflektierte Strömung anströmen müßte. Die Strömung wird aber in 
ihrem Überschall- oder in ihrem Unterschallteil physikalisch bis nahe an den Grenzkreis 
realisierbar, wenn man für einen geeigneten räumlichen Zu- bzw. Abfluß sorgt. Andernfalls 
wird die Strömung ihren adiabatischen Charakter, sei es stetig oder stoßweise (Verdichtungs- 
stoß) schon vor dem Grenzkreis verlassen oder eventuell rückläufig ihren Anfangszustand ver- 
ändern. Die anderen Beispiele dieser Arbeit werden weitere Einsicht in die Natur des Strö- 
mungsstoßes gewähren 

Zunächst entsteht die Frage, ob es unter der Voraussetzung, daß die Strömungs- 
geschwindigkeit bis zur Maximalgeschwindigkeit wächst, überhaupt Strömungen gibt, die frei 
von Strömungsstößen sind, abgesehen von der Parallelströmung. Wir gelangen zu einer ge- 
suchten Strömung, wenn wir nach allen Strömungspotentialen Y fragen, die nur 
von w abhängen. Unsere Gl. (15) liefert dann für Y die Bedingung 


wobei 
w 
gesetzt ist. 
Die Lösung dieser Gleichung lautet 
dw 
w 
oder ausintegriert 
1 3 5 
Ii+l-am) 


Es folgt daher mit Benutzung von (13) 


*%) Im Sinne der in der Einleitung angegebenen Definition. 
®) Eine allgemeine Theorie des Strömungsstoßes erscheint an anderer Stelle. Es sei hier nur gesagt, daß die 
Stoßlinien sich allgemein als natürliche Grenzen einer adiabatischen Strömung erweisen. Sie genügen der Bedingung 


1 1 2 


| 
ds 
-5 
Auf dem Grenzkreis ist aber w=e, also 
| 
dw_dw 1 
ds dt ds w’ 
dt 
5 
1— am?) 
an’). 
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Die Gl. (11) ergeben 


sin d cos d 
sin d 
de 
Hieraus findet man für x und y selbst, von additiven Integrationskonstanten abgesehen, 
sin 
w w 


Für die Stromlinien gilt wieder Y/ = const, also w=const=k. Die Stromlinien lauten 


die zugleich die Kurven konstanter Geschwindigkeit sind. 


Es handelt sich hier um einen Wirbel'’), der die gleiche Geschwindigkeitsverteilung 
aufweist wie der inkompressible Wirbel, nur daß das Innere des Wirbels von einer Kreis- 
fläche gebildet wird, auf deren Rand die Maximalgeschwindigkeit und damit die Dichte Null 
erreicht wird. Ein Strömungsstoß tritt hier also nicht auf. 


Nach den Ausführungen von S. 188 erhalten wir aus der Quell- und der Wirbelströmung 
sofort eine neue Strömung, die wir in Analogie zur inkompressiblen Strömung als Wirbel- 
quelle (bzw. -senke) bezeichnen. Der Zusammenhang zwischen x, y und w, ® lautet 


sin d 0, 
0, sin d 
— — — 
1 m 29 w 


Für die Stromlinien gilt die Beziehung 


ı 
„2)2 1 
“- 1+l1-am?) 


- a +0,9—k 


(k konstant), wobei wieder 


zu setzen ist. 
Dieses Beispiel (Bild 2)'*) 
läßt erstmals den allge- 
meinen Charakter 
einer adiabatischen 
Gasströmung erken- 
nen. Man beobachtet näm- 
lich auf jeder Stromlinie 
einen als Spitze in 
Erscheinung treten- 
den Strömungsstoß. 


Bild 2. 


3. Strömung um eine Kante. Die Differentialgleichung (15) 


mit 
„—1 +1 
10, G. J. Taylor: Recent work on the flow of eompressible fluids, Journal of London Mathem. Soc. Bd. 5 (1930), 
S. 224; Deutsche Übersetzung: Z. angew. Math. Mech. Bd. 10 (1930), S. 334. 
11) Bei der Durchführung der Rechuungen und Zeichnungen zu dieser Arbeit wurde der Verfasser in daukens- 
wertester Weise von den Herren A. Fendt und H. Klee unterstützt, 
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soll nun allgemeiner integriert werden. Wir suchen zunächst Lösungen von der Form 


wo P eine Funktion von w und Q eine Funktion von ö allein sein soll. Haben wir beliebig 
viele Lösungen dieser Art, so ergeben ihre linearen Kombinationen mit konstanten Koeffizienten 
weitere Lösungen von sehr allgemeiner Art. Die Ableitungen von P und @ nach ihren 
Argumenten werden mit Strichen bezeichnet. Dann wird auf Grund der Differentialgleichung 


Pin) Pin) ' 00) 
Da die ersten beiden Summanden lediglich von w abhängen, der dritte aber nur eine Funktion 
von ® ist, müssen die Beziehungen 


+ =0. 


gelten, wo A eine Konstante ist. Wir beschränken uns auf ein positives 
Dann ergibt die Integration der zweiten Gleichung 
(C, ®, Integrationskonstanten). Die erste Gleichung lautet 
.... (Al). 


Wir betrachten nun den besonders einfachen Falln—1: Ein partikuläres Integral der 
Differentialgleichung ist jetzt 


Po)=— 
die allgemeine Lösung daher 
1 1-?w?dw 
(C,, €, Integrationskonstante). 


Von besonderem Interesse ist die spezielle Lösung der Differentialgleichung (15): 


sin ®, 
w 
Wie in den übrigen Beispielen berechnen sich hierfür die Werte 
C 
sind, 0050, 
2 
cos’ 0, sin2# 
sin29 sin®cosd 0, 60829 
Hieraus folgt 


Der Wert des Integrals ist 


1 [1 _ 1 1 1 1 1 
am’)? (1 — am”)? (1— am?)?® 
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Für die Stromlinie gilt wegen Y=const: 


w—ksin®. 


Führt man dementsprechend in den Gleichungen für x und y an Stelle von # das w ein, so 
findet man für die Stromlinien die Parameterdarstellung 


y=+c® 


Dabei ist J die durch (24) gegebene Funktion von w. 


Der Charakter der Strömung ergibt sich sofort, wenn man sie sich zunächst inkom- 
pressibel denkt. Setzt man dementsprechend e=c,=x, o=o,. so vereinfachen sich die 
Gl. (25) zu 


— 
Dabei haben wir über die Integrationskonstanten so verfügt: 
C=1l, 


Das bedeutet keine wesentliche Beschränkung der Allgemeinheit. Eliminiert man aus den 
Gl. (26) noch das w, so wird 


1 


Diese Gleichung stellt, wenn 
man k verschiedene konstante 
Werte erteilt, eine Schar von 
konfokalen Parabeln dar, 
in denen die Strömung um 
eine Halbgerade (Kante) 
erfolgt (Bild 3). 

Die kompressible Strö- 
mung hat für kleine Geschwin- 
digkeiten nahezu die gleiche 
Form (vgl. Bild 4). Z. B. ist die 
Stromlinie k=300 der Zeich- 
nung noch von parabolischer 
Gestalt. Die höchste Geschwin- 
digkeit wird auf der Symmeetrie- 
geraden angenommen. Sie stimmt 
mit dem k-Wert überein. Auf 
der Stromlinie k = 300 ist also 
die Geschwindigkeit im Scheitel 
gleich 300. Mit Annäherung an 
die Kante wird die Geschwin- 
digkeit höher, die Kurven wei- 
chen von der Parabel mehr und 
mehr ab, indem sie in der Um- 
gebung des Scheitels flacher 
werden. 

Die Kurven, auf denen die 
Geschwindigkeit konstant ist, 
sind die Kreise 


Bild 3.%.Strömung um eine Kante (inkompressibel). 


Sie sind in Bild 4 eingezeichnet 
und durch die auf ihnen ange- 
nommenen Geschwindigkeiten 
gekennzeichnet. y Bild 4. Strömung um eine Kante (kompressibel). 


XUM 


N 
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Zur Diskussion der Stromlinien berechnen wir noch die Tangente und die Krümmung 4 . 


Man findet aus den Kurvengleichungen (25) 


da k ( 

und 

1 9a k 


Auf der Kurve k ist die Geschwindigkeit höchstens gleich k und natürlich höchstens gleich 


der Maximalgeschwindigkeit. Für w=k wird =». Im übrigen ist SU durchaus stetig. 


Ist k größer als die Maximalgeschwindigkeit, so ist I: überall stetig. Von besonderem Inter- 
esse ist nun der Nenner von =. Ist dieser nämlich Null,. ist also 


w? 
so hat die Kurve trotz stetiger Tangente unendliche Krümmung, also eine Spitze. Berechnet 
man mit Hilfe der Formel 
b dw 


die Beschleunigung 5, so ergibt ‚sich 


(32). 


Der Nenner ist derselbe wie bei . wird also Null. Der Zähler nimmt wegen (30) den Wert 


° ew==0 an. Somit ist an der Spitze die Beschleunigung 5b unendlich, d. h. 
es liegt ein Strömungsstoß vor. 
Um auf einer Stromlinie k den evtl. vorhandenen Stoß zu ermitteln, drücken wir in Gl. (30) 
c® durch w? aus und erhalten zur Bestimmung des zum Stoß gehörigen w eine quadratische 
Gleichung für w?, aus der sich 


EIN 
ergibt. Ein Stoß ist demnach dann und nur dann vorhanden, wenn 
2 
also 
34) 
für k=1,4 somit 


ist. 

Da k zugleich die Höchstgeschwindigkeit auf einer Stromlinie angibt, haben wir das 
folgende für die Praxis wichtige Ergebnis: Bei einer Strömung über eine gewölbte 
Fläche ist die Schallgeschwindigkeit im allgemeinen kein kritischer Zu- 
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stand. Kritisch sind nur Strömungsstöße, die im allgemeinen nicht mit, 
sondern erst oberhalb der Schallgeschwindigkeit auftreten. Im vorliegenden 
Beispiel tritt der Stoß auf der Stromlinie k = 1,667 c, etwa bei der Geschwindigkeit w=1,2 c, 
auf. Eliminiert man aus (25) und (30) den Parameter k, so erhält man den Ort") aller 
Strömungsstöße (die Stoßlinie) in der Form 


J 


mw’ 


Insgesamt zeigt unser Beispiel (Bild 4) folgendes: Bei kleineren Ge- 
schwindigkeiten sind die Stromlinien näherungsweise konfokale Parabeln. 
Die Stromlinien zeigen auch bei größeren Geschwindigkeiten einen glatten 
Verlauf, wenn nur die auf der Stromlinie erreichte Höchstgeschwindigkeit 
kleiner als 1,667c, ist. Wird die Geschwindigkeit 1,667c, erreicht, so 
treten an den Stellen mit etwa 1,2c, Geschwindigkeit erstmals Strömungs- 
stöße auf. Auf jeder Stromlinie, auf der eine Geschwindigkeit >1,667c, er- 
reicht wird, liegen 4 Strömungsstöße, welche Spitzen der Stromlinie sind. 
Physikalisch sinnvoll ist immer nur der Teil der adiabatischen Strömung, 
der keine Strömungsstöße überschreitet. Z.B. also istder Teilder Strömung 
sinnvoll, der durch die ausgezogenen Stromlinien dargestellt ist. Die Ge- 
schwindigkeit der Strömung im Unendlichen ist Null. Bei Eintritt der ört- 
lichen Schallgeschwindigkeit ist die Entfernung der Stromlinien nahezu 
ein Minimum. Die Stellen mit örtlicher Schallgescehwindigkeit sind dureh 
Punkte auf den Stromlinien gekennzeichnet. Die Geschwindigkeitskom- 
ponente senkrecht zur Stoßlinie ist gleich der lokalen Schallgeschwindig- 
keit (folgt aus (30) und (36)). 


4. Allgemeines Integrationsverfahren. Anwendung auf die Doppelquelle. Im vorigen Para- 
graphen haben wir die Differentialgleichung (21) 


P"+w(l Bw) ...... 87), 


auf deren Lösung es ankommt, für den Fall n=1 integriert. Nunmehr soll dureh 
Reihenentwicklung eine Lösung P, für beliebiges n >0 gefunden werden’). 
Die Theorie der linearen Differentialgleichungen zeigt, daß im vorliegenden Falle der 
Ansatz 


zum Ziele führt. Differenziert man zweimal nach w und setzt man die gefundenen Reihen 
in die Ausgangsgleichung ein, so ergeben sich durch Koeffizientenvergleich zur Bestimmung 
der unbekannten ce; folgende Gleichungen: 


und wenn man zur Abkürzung 
setzt, 
(beliebig), 


Die P,„ konvergieren bis zur Maximalgeschwindigkeit (ausschließlich). 
Für Y erhält man nach den Ausführungen von S. 192 z. B. die Lösungen 


12) Ju Bild 4 ist nur ein Teil dieses Ortes klein gestrichelt eingezeichnet, der für den praktischen Verlauf der 
Strömung von Interesse ist. 
13, Ein ähnliches aber wesentlich komplizierteres Verfahren wendet C. Ferraria.a. O0. an. 


Allgemein gilt die sehr einfache Rekursionsformel: 
XUM 
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Der Zusammenhang zwischen x,y und w,d ergibt sich ebenso wie in den früheren Spezial- 
fällen. Zunächst hat man für die erste Lösung I 


w PR = = gilt. Die Beziehungen (13) liefern 


Die Gl. (11) ergeben dann nach einiger Rechnung für x und y die folgenden mit »„ und 9, 
bezeichneten Werte: 


) 1 „tn cos(n—1) 2» 
bzw. für die zweite Lösung Y (46) 
P, 
(n=F1), 
und im Falle n=1 
1 
wobei 
1(o 
am, 
bzw. 
(Pr + + K,, 
wobei 


gilt (K,, K, Integrationskonstante). Die Stromlinien sind gegeben durch 


(k konstant). Nach dem Überlagerungsprinzip (S. 188) lassen sich diese Lösungen linear zu 
neuen Lösungen kombinieren. Bildet man auf diese Weise unendliche Reihen, so müssen 
diese natürlich gleichmäßig konvergieren. 


Gehen wir zu inkompressiblen Strömungen über, so wird 


Will man also eine kompressible Verallgemeinerung irgendeiner inkom- 
pressiblen Strömung herleiten, so hat man für diese das Strömungs- 
potential Yin der Form _ 


B„sinnd) 


zu entwickeln und alsdann w” durch P,„ zu ersetzen. 
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Ein Beispiel möge dieses Verfahren erläutern. Wir suchen eine kompressible Ver- 
allgemeinerung der Doppelquelle. Setzt man z=x+iy=reir, so lautet die Strömungs- 
funktion für die inkompressible Doppelquelle im Nullpunkt 


abgesehen von einer unwesentlichen multiplikativen Konstanten. Die Geschwindigkeits- 
funktion lautet 


Es gilt also 
1 
und daher 
1 


Y ist der Imaginärteil von {. Man hat somit 


Für die kompressible Verallgemeinerung ist 


v__ Pı cos 


zu setzen, wobei für Pı die Entwicklung 


1 
Pı=w: lei 
2 2 


2 ... 


gilt. Nach den Formeln (39) bis (42) findet man, wenn man Ü=1 setzt, 


—= 0,000 173 & 


Wir beschränken uns auf drei Glieder und haben 


(1 (2) + 0,00969 ()) 


8 0 
1 9 
Pr 0000 


Die Strömung ist somit gegeben durch 


1 
=k. 
2 
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Diese Strömung ist in Bild 6, die zugehörige inkompressible in Bild 5 gezeichnet. Auch dieses 
Beispiel zeigt grundsätzlich das gleiche Verhalten wie die früheren Beispiele. 


Bild 5 (links). 
Doppelquelle 
(inkompressibel). 


Bild 6 (rechts). 
Doppelquelle 
(kompressibel). 


Die obigen allgemeinen Betrachtungen mögen schließlich noch durch einige Bemerkungen 
über die Machschen Linien ergänzt werden. Es sind das die Kurven, senkrecht zu denen 
die Geschwindigkeitskomponente der Strömung gleich der lokalen Schallgeschwindigkeit ist. 
Ist r der Neigungswinkel der Tangente einer Machschen Linie, so gilt allgemein 


sin 


Hieraus folgt 
2 2 
(1+(2) eos — sin 


w? dx dx 
und mit Benutzung von (11) 


mw 


Die Lösung dieser Differentialgleichung lautet 


Das in diesem Paragraphen dargelegte Integrationsverfahren ist olıne weiteres zur Be- 
handlung von Profilströmungen geeignet. Von besonderem Interesse wird dabei die 
Untersuchung der Lage der Stoßlinie sein und die Klärung der Frage, wie ein Profil gestaltet 
sein muß, damit ein Strömungsstoß trotz möglichst hoher Geschwindigkeit nicht auftritt. 
Diesen Fragen ist eine Fortsetzung der vorliegenden Arbeit gewidmet. 


=+ 


Der Verfasser möchte schließlich Herrn A. Lippisch für manche Anregung herzlich 
danken, die er durch gemeinsame Diskussionen der hier behandelten Probleme erhalten hat. 
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Über das Kerbproblem in der Plattentheorie. 
Von H. Neuber in Braunschweig. 


1. Einführung. Die bisher vorliegenden Lösungen des Kerbproblems') haben die Möglich- 
keit gegeben, die Verteilung der Spannungen in verschiedenartig geformten Bauteilen rech- 
nerisch zu erfassen. Es fehlt jedoch noch bezüglich plattenartiger Bauglieder an Rechnungs- 
grundlagen. Um Mißverständnissen vorzubeugen, sei hier ausdrücklich erwähnt, daß mit der 
Bezeichnung „Platte“ ein flacher Körper mit ebener Mittelfläche gemeint ist, der — im Gegen- 
satz zur „Scheibe* — nur durch solche Kräfte belastet ist, die senkrecht zur Mittelebene 
wirken, bzw. durch Momente, deren Drehachsen der Mittelebene parallel laufen. Vorliegende 
Arbeit bringt strenge Grundlösungen des Problems der gekerbten Platte und soll dazu bei- 
tragen, die in der Kerbspannungstheorie bezüglich der Platten noch bestehende Lücke aus- 
zufüllen. Im Ergebnis werden Beziehungen zwischen Randform und Spannungsverteilung 
gewonnen, welche zur Beurteilung der Kerbwirkung bei Platten aufschlußreiche Grundlagen 
liefern. 


Il. Ableitung einer Plattentheorie aus dem Dreifunktionenansatz. Die meist angewandte 
Plattentheorie ist die Kirchhoffsche?); sie geht von der Annahme aus, daß die Biege- 
spannungen über die Plattendicke linear verteilt sind. Eine genauere Formulierung, welche 
den elastischen Grundgleichungen genügt, gab Michell?); hiernach sind die Biegespannungen 
nicht proportional dem Abstand z von der Plattenmittelfläche, sondern enthalten außer z noch 
2°. Bei Aufstellung der Ausdrücke für Biegemomente und Schermoment durch Integration 
der mit z multiplizierten Spannungen o,, o, und r,, über die Plattendicke tritt außer dem 
Faktor ö° noch Ö° auf. Das Auftreten von 2° bei den Spannungen und ö° bei den Momenten 
bringt in zweifacher Hinsicht eine Abweichung von der Kirchhoffschen Theorie mit sich, 
und es ist in der Michellschen Darstellung nicht klar zu übersehen, von welcher Größen- 
ordnung diese Abweichungen sind, da das Glied mit ö° bei den Momenten mit in die Rand- 
bedingungen eingeht. Es erscheint mir deshalb angebracht, an Stelle der Michellschen 
Lösung eine andere aufzustellen, welche den Vorzug besitzt, daß bei den Momenten Ö° 
nicht mehr auftritt. Das Auftreten von ö° in der Michellschen Lösung liegt in der 
Verwendung der Durchbiegung der Plattenmittelfläche als Spannungsfunktion begründet. 
Wie ich nachstehend. zeigen werde, läßt sich mit Hilfe des Dreifunktionenansatzes*) eine den 
elastischen Grundgleichungen entsprechende Formulierung derart angeben, daß bei den 
Momenten nur noch der Faktor ö° allein auftritt. Diese Darstellungsart hat den Vorteil, 
daß bei den Spannungen die Anteile mit z unmittelbar den Momenten entsprechen, während 
die zusätzlichen Anteile mit 2° und ö°z für sich innerhalb der Seitenfläche eines Platten- 
elementes ein Gleichgewichtssystem bilden, dessen Einfluß auf die Spannungsverteilung der 
Umgebung nach St. Venant nicht wesentlich sein kann. Die so vereinfachte Theorie der 
dicken Platte stellt sich mittels des Dreifunktionenansatzes wie folgt dar. 

Es sei F,(x,y) eine nur von den Koordinaten x, y der Plattenmittelebene abhängige 
biharmonische Funktion, die sich aus zwei ebenen harmonischen Funktionen y, und y, in 
der Form 


zusammensetzen möge. Dann ist 


eine harmonische Funktion. Weiter überzeugt man sich leicht, daß 


2 
F—Z4F, und — AR, 


räumliche harmonische Funktionen sind, also der Potentialgleichung d ®=0 genügen. Man 
erhält nun mit Hilfe dieser Funktionen in einfacher Weise die geforderte allgemeine Lösung 
für die dieke Platte, wenn über die vier harmonischen Funktionen ®,, ®,, ®, ®,, welche in 
meinem für alle Probleme der Elastostatik geltenden Dreifunktionenansatz*) auftreten, 


1) H. Neuber: Kerbspannungslehre, Berlin 1937, Verlag Julius Springer. 
2) G. Kirchhoff: Crelles Journ., Bd. 40 (1850), S. 51; siehe auch A. u. L. Föpp], Drang und Zwang, Bd. 1, 
München und Berlin 1924, $ 17. 

3) J. H. Michell: Proc. Math. Soc. London, Bd. 31 (1899), 8. 100; ferner A. E. H. Love: A treatise ‘on the 
mathematical theory of elastieity, deutsch von A. Timpe, Leipzig 1907, $$ 303 u. 304. 
4) Siehe 1); ferner H. Neuber: Z, angew. Math. Mech. Bd. 14 (1934), S. 203 bis 212. 
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folgendermaßen verfügt wird: 


1 3°? 1 1 1 


D,—0, d,—0, 


1 


Die Konstante a hängt mit der Poissonschen zusammen und berechnet sich aus 


| 


m 


Entsprechend der Theorie des Dreifunktionenansatzes*) bestimmt man nun zunächst eine 
biharmonische Funktion F aus der Gleichung 


und erhält die Komponenten des elastischen Verschiebungsvektors aus 

D,, 

36. = a®, 


Für die Spannungen o,, 7x, usw. gelten andererseits die Beziehungen 


(7). 
od, 
Hier ergibt sich durch Einsetzen von (3) in (5) 
Für die Verschiebungskomponenten folgt 
1 9 
Schließlich ergeben sich folgende Ausdrücke für die Spannungen: 
ar) 
1 
1 
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Man erkennt, daß es sich in der Tat um den Spannungszustand einer Platte handelt: Beider- 
seits lastfreie Oberfläche; denn für 2= + 3 verschwinden die Schubspannungen r,, und ryz, 


ferner tritt o, nicht auf (bei belasteter Oberfläche kämen noch weitere Glieder hinzu). Die 
Spannungen o, und o, liefern die Biegemomente, die Schubspannung r,, das Schermoment, 
schließlich die Schubsspannungen r,, und r,, die Querkräfte. Weiterhin erkennt man, daß 
die Glieder mit 2° und ö?z von der Michellschen Lösung abweichen; sie treten hier in 


der Form 
6? 2° 


auf. Man erkennt, daß die zur Bildung der Momente erforderliche Integration der mit z 
multiplizierten Spannungen o,, o, und r,, für diese Anteile keinen Beitrag liefert. Die 
Glieder höherer Ordnung beeinflussen also bei der hier gegebenen Dar- 
stellung nicht den Kraftfluß, damit auch nicht die Randbedingungen. Mit- 
hin gehen sie auch nicht in die für die Konstantenermittlung maßgeblichen Gleichungen ein. 


III. Übergang auf krummlinige Koordinaten. Zur Behandlung von Kerbproblemen ist die 
Verwendung krummliniger Koordinaten erforderlich; nur bei Durchführung der Rechnung mit 
einem Bezugssystem, welches dem Rand angepaßt ist, besteht die Möglichkeit, den Krüm- 
mungsparameter allgemein mitzuführen. Das Bezugssystem sei festgelegt durch die Gleichungen 


Den Koordinaten u,v,w mögen die Indizes 1, 2, 3 entsprechen, dann erhält man die neuen 


Spannungskomponenten entsprechend den für den Dreifunktionenansatz in gekrümmten Ko- 
ordinaten geltenden Beziehungen®) in folgender Form: 


-Zurl, 


F* 


2 
- 


Hierin wurde zur Abkürzung 


a 


gesetzt. Für die in (12) auftretenden Operatoren gilt 


14 


wobei 


5) Vgl. 


_ 
h,dulhdu) 
h,öv) Mn,dudu’ 
h,dulh,dv/) hih,dvdu 
| 
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IV. Vereinfachte Darstellung der Randbedingungen. Für die Durchführung der Rechnung 
bei strenger Einhaltung der Randbedingungen bedeutet es eine wesentliche Erleichterung, 
wenn sich eine Vereinfachung jener Gleichungen vornehmen läßt, welche die Randbedingungen 
darstellen. So läßt sich z. B. der ebene Spannungszustand der Scheibe in einfacher Weise 
als Randwertaufgabe der Airyschen Spannungsfunktion auffassen, derart, daß in den Rand- 
bedingungen nur mehr die ersten Ableitungen der Spannungsfunktion auftreten. Diese Ver- 
einfachung beruht auf einer Integration der betreffenden Gleichungen längs des Randes. Wie 
ich nachstehend zeigen werde, läßt sich eine ähnliche Vereinfachung auch in der Platten- 
theorie durchführen. 

Bild 1 zeigt eine Plattenecke, 
welche längs der neuen Koordinaten- 
flächen herausgeschnitten ist. Die rechte 
Schnittfläche (v=konst.) sei zugleich 
Begrenzung der Platte. Hier muß, da 
sich die Randbedingungen bei Kerb- 
problemen in der Regel auf die Span- 
nungen und nicht auf die Form- 
änderungsgrößen beziehen, das Biege- 


moment 
d 
T 
/ M,=lo,sds. ... . (ie 
Resthraft:/t,2dı Jtydzh,du 
-f f 2 
Bild I. Kräftespiel an einer Plattenecke. einen vorgegebenen Verlauf annehmen. 


Andererseits haben auch die übrigen, 
am Plattenrand angreifenden Spannungen eine Bedingung zu erfüllen. Man erhält 
diese zweite Bedingung, indem man aus dem längs Ah, dw wirkenden Schermoment 


\7,2dz h,du jeweils ein Kräftepaar bildet, wobei die zugehörigen Ersatzkräfte in der z- 


Richtung wirken (Bild 1)°). Durch Zusammenfassung dieser Ersatzkräfte für benachbarte 


Elemente mit der Querkraft |r,,dzh,du ergibt sich schließlich eine auf die Längeneinheit 


des Randes bezogene resultierende (Querkraft vom Betrage 


2 2. 


Bei Anwendung von (12) und (14) gehen (16) und (17) über in 


ö° 2 

2 
Am lastfreien Rand gilt demnach 

2 
5 

2 


Die nun folgende Umformung dieser Bedingungen beruht auf geometrischen Zusammen- 
hängen zwischen den Ableitungen nach x und y einerseits und den Operatoren [ ],, [ ],, usw. 
andererseits, welche sich daraus ergeben, daß sich die ersten Ableitungen wie Vektor- 


6) Die Zulässigkeit einer solehen Umwandlung von Torsionsschubspannungen in zusätzliche Querkräfte am 
Plattenrand wurde erstmalig von Thomsen und Tait nachgewiesen; vgl. Handb. d. theor. Phys., Bd.1, T.2, Art. 645 ff. 
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komponenten, die zweiten wie Tensorkomponenten transformieren. Wird der jeweilige Winkel 
zwischen den Richtungen & und 1 bzw. y und 2 mit @ bezeichnet, so gilt: 


11 12 (4) 
Ferner bestehen zwischen AF, 23H * und der konjugierten Funktion — am die Cauchy- 
Riemannschen Gleichungen 
so daß (21) übergeht in 
_ 
oder nach Integration längs des Randes 
4 od y, 
Andererseits folgt aus (20) 
4oy, 
(28). 
Setzt man diese Ausdrücke in (24) ein, so ergibt sich 
F, N) 
oder wegen (22) und (23) 
4 
4 


Hierbei ist die neue Funktion y, die zu y, konjugierte Funktion, erfüllt also die Gleichungen 


dy 


Die Beziehungen (31) und (32) lassen sich, wie man sieht, ebenfalls längs des Randes 
integrieren. Dann ergibt sich 


(33). 


ar, 4 
d 4 


Hieraus folgen 


oF, 4 dx dx 


14* 
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Diese einfachen Beziehungen, welche in derselben Form auch für einen Rand «= konst. 
gelten, ermöglichen die Anpassung der Lösung an die vorgeschriebenen Randbedingungen, 
ohne zunächst die Aufstellung der meist verwickelteren Ausdrücke für die Spannungen er- 
forderlich zu machen. Dadurch vereinfacht sich die Ermittlung der im Lösungsansatz ent- 
haltenen Konstanten ganz wesentlich, wie insbesondere aus den nachstehend behandelten 
Kerbproblemen hervorgeht. Zuvor sei auch für den Kraftfluß eine vereinfachte Darstellungs- 
form angegeben. 


V. Auflösung der Kraftflußintegrale. Es handle sich um den Kraftfluß, der durch eine 
Schnittfläche « = konst. hindurchgeht, welche zwei lastfreie Ränder (Parameter v=v, und 
®=v,) miteinander verbindet. Unter Berücksichtigung der an den lastfreien Rändern noch 
übrig gebliebenen Restkräfte (Bild 1) wird das um die x-Achse drehende Moment der 
Schnittkräfte 

2 


v2 


2 


Mit Hilfe von (24) und (25) läßt sich der Ausdruck auf folgende Form bringen: 


oder 


cosp dv . (39). 


410 1 4 1 
oder mittels (23) 
v2 
ö° 4 


Entsprechend (27) und (34) gilt für den lastfreien Rand v=»;: 


4dy\ 4 
Ebenso, jedoch mit einer anderen Konstante, gilt für den lastfreien Rand v=»,: 

OF, 

Demnach folgt aus (41) die einfache Beziehung 

6° 
(44). 


In derselben Weise läßt sich nachweisen, daß folgende einfache Beziehung für das um die 
y-Achse drehende Moment gilt: 


M=&8B,-B). m. 


Schließlich ist noch die im Schnitt übertragene Gesamtquerkraft anzugeben, die mit P; 
bezeichnet sei. Bei Berücksichtigung der Restkräfte ergibt sich 


d d 
v2 
| t,h,dvdz+ | — 2dz 


Nun gilt gemäß (27) am Rand v=»;: 
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Am Rand v=v, besteht dieselbe Beziehung, jedoch mit einer anderen Konstanten, die mit C, 
bezeichnet sei. Dann folgt aus (46): 


6? 


Durch die Gl. (36), (37), (44), (45) und (48) sind die Randbedingungen und der Kraftfluß der 
Platte in einfacher Form festgelegt, und ich glaube, hiermit für das Randwertproblem der 
Platte allgemein eine beträchtliche Vereinfachung erzielt zu haben, zumal es nunmehr — wie 
beim ebenen Spannungszustand — möglich ist, die Spannungsfunktion zu bestimmen, ohne 
zunächst die meist unübersichtlichen Ausdrücke für die Spannungen aufstellen zu müssen. 
Dieser Vorteil tritt bei den nachstehend behandelten Kerbproblemen besonders hervor. 


VI. Die Grundprobleme der gekerbten Platte. Uın für die Berechnung der Spannungen 
in gekerbten Platten eine Grundlage zu schaffen, ist es erforderlich, jene Grundlösungen auf- 
zustellen, die auch bei der Scheibe (Zug- oder Biegestab) den Ausgangspunkt der Kerbspannungs- 
theorie bildeten: Beiderseitige Außenkerbe (A) und Langloch (B). Die Lösungen 
wurden durch Idealisierung zur hyperbolischen bzw. elliptischen Randform gewonnen, und 
ihre Anwendung auf praktische Fälle mit beliebigen Abmessungen durch Interpolationsformeln 
sichergestellt”), die auf der Grundlage des Abklingungsgesetzes entwickelt wurden. Nachdem 
die Gültigkeit dieser Interpolationsformeln in gleicher Weise auch für die Platte vorausgesetzt 
werden darf, handelt es sich ausschließlich um die Aufstellung der strengen Lösungen dieser 
Grundprobleme für die Platte (Bild 2 und 5). Als Bezugssystem dienen die elliptischen Ko- 
ordinaten mit 


z=©inucosv, y=6ojusinv, 


(49). 


h,=h,=h=y&in’u-+ cos’ v = (Cof 2u+ cos2v) 


Bei Berechnung der Spannungen sind in den folgenden Abschnitten die für dicke Platten in 
Betracht kommenden höheren Glieder noch unberücksichtigt geblieben (d.h. es wurde mit F, 
statt F* gerechnet), um möglichst einfache Endformen zu gewinnen. Die Durchführung der 
zugehörigen Ergänzungsrechnung ist jedoch im Rahmen der Behandlung weiterer Platten- 
probleme zu einem späteren Zeitpunkt beabsichtigt. 


A. Erstes Grundproblem: Platte mit beiderseitiger tiefer Außenkerbe 
bei Biegung. Gemäß Bild 2 handelt es 


sich um eine biegebeanspruchte Platte mit ı 2? 
tiefer beiderseitiger Außenkerbe. Der Rand C BR 2) 
M 2 M 


sei durch die Hyperbel ve=+v»,,d.h.v,=v,, 
v‚=—v, festgelegt. Die Randbedingungen 
(36) und (37) haben für das gewählte System 4 


die Form (mit ve=: v,): Bild 2. Gebogene Platte mit beiderseitiger Außenkerbe. 


(— y, &oj u cos v, — y, Sin u sin v,) 


— (,sinv,cosv, + A, cosv,—+ B, ©&in u sin v,, 


(50), 
Eof u cosv,—+ y, Sin u sin v,) 
= (, sin v, A, Eof u cos v, — B, Sinusinv, 
+ (y, Sin u sin v, — y, Cof u cos v,) 
—= (,&inu Cofu + A,Sinusinv, + B, Cofu cos v,, (51) 


+ (—y, Sin u sin v, — y, Cof u cos v,) 


—=(, Sin u Eof u — A, Sinusinv, + B, Cofucosv, 


Von den für den Kraftfluß maßgeblichen Größen M,, M,, P; ist nur M, von Null verschieden. 
Entsprechend (44), (45) und (48) wird 


?) Kerbsp gslehre, S. 5. 
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Die Lösung für F, muß gemäß 

aus zwei harmonischen Funktionen y, und y, hergestellt werden, wobei obige Bedingungen 
erfüllt sein müssen. Außer y, und y, tritt auch die zu y, konjugierte Funktion y, auf. Als 
geeignete Lösung kommt allein die folgende in Betracht: 


y,— A Cof u cos» + B (v Eof u sin — u Sin u cos v), 
y‚=Cu, y‚=Cv 
Es wird hiermit 
BJu&inucose . ..... (BD). 


Die Erfüllung der Randbedingungen verlangt, wie man unmittelbar erkennt, daß die Kon- 
stanten A,, A,, C,, €, verschwinden und BB=—B,=6M/ö* zu setzen ist. Gl. (50) ver- 
langt ferner 


Sinu(Acosv, + Bv, sin v,) + (C — B) (u Eoj u+ Sinu)cosv, 


C (— u u cosv, — v, Ein u sin v,)= Sinusinv®, 
Hieraus folgen 
Andererseits folgt aus (51) 
u (— Asinv, + B(v, cos v, + sin v,)) — (C— B) u Sinusinv, 
+ (u Sin u sin v, — v, Eof u cos v,) = Cojucosv, we. 
Demnach müssen folgende Gleichungen erfüllt sein: 
— Asinv, + B(v, cosv, + sin — Cv,cosv,—= 


Gl. (60) ist mit (57) identisch. Die Auflösung der drei Gl. (57), (58) und (61) ergibt für die 
drei Konstanten A, B, C folgende Werte: 


B=c(1-), 
a a 
121 


Durch Einführung der elementaren Biegespannung 


3M 
der 
wird wegen a=sinv,: 


ö sine». 


Die Spannungskomponenten ergeben sich aus (12) bis (15). Bei Berücksichtigung der für die 
Konstanten gewonnenen Werte erhält man: 


Z. angew. Math. Mech. 
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Die Glieder mit z ö° und 2° bei o,, o, und r,,, die den Kraftfluß nicht beeinflussen (s. Abschn. II!), 
sind hierbei unterdrückt. 

Man überzeugt sich leicht, daß die gewonnenen Spannungskomponenten in der Tat die 
Randbedingungen erfüllen. Die Spannungsverteilung zeigt die für Kerbprobleme typischen 
Erscheinungen (starke Spannungsspitze im Kerbgrund, die mit der Krümmung der Kerbe 
zunimmt, Spannungsabfall nach dem Inneren und längs des Randes). Bild 3 erläutert den Verlauf 
der Spannungen im Symmetrieschnitt und längs des Randes für tge,=y3. Mit Einführung 
des Krümmungshalbmessers o ergibt sich für die Höchstspannung folgender Ausdruck: 


l\ı/a 
> 
Imax__ V (66) 


[3 + arctg V — “zZ 


Die Abhängigkeit von V - ist in Bild 4 ersichtlich. 


5 


PH Rh: 03; ne Bild 4. Die Höchstspannung in gebogenen Platten mit beider- 
seitiger Außenkerbe in Abhängigkeit von der Kerl'krümmung. 


L Bild 3. Bild 3. Spannungsverlanf in einer gebogenen Platte mit 
beiderseitiger Außenkerbe. 


Vergleicht man das Ergebnis mit dem entsprechenden Scheibenproblem*), so zeigt sich, 


T — 0,394 (tür 0,3) multi- 
plizierten Formzahl der Scheibe ist. Bei dem nachstehend behandelten zweiten Grundproblem 
werde ich zeigen, daß derselbe Zusammenhang besteht. 

B. Zweites Grundproblem: Die Platte mit Langloch bei Biegung. Ent- 
sprechend Bild 5 handelt es sich jetzt um eine gebogene Platte mit Langloch, wobei ellip- 
er tische Lochform vorausgesetzt wird. Ferner 

ist angenommen, daß der Abstand von den 
7 äußeren Begrenzungen der Platte genügend 


daß die Formzahl der Platte im wesentlichen gleich der mit 


groß ist, so daß die Kerbspannungen am 

Außenrand nicht mehr in Erscheinung 

e treten. Die Lösung enthält zunächst jene 

177 Glieder, welche einer ungelochten Platte 
Bild 5. Gebogene Platte mit Langloch. bei reiner Biegung entsprechen: 

s) Kerbspannungslehre, S. 36. 


XUM 


2 9\Eojusine 2cos’r 
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— 2a) (y’ — x), 


pa 
2p 


Der Spannungsbezugswert p entspricht der Biegespannung in der äußersten Faser. Die Lösung 
für die gelochte Platte wird erhalten, indem dieser Ansatz in den elliptischen Koordinaten 
geschrieben wird und geeignete Zusatzfunktionen herangezogen werden, welche für große Ent- 
fernung vom Loch verschwindend kleine Werte liefern. Auf diese Weise erhält man: 


Sinu —4Be*) cosv, 
(68). 


2a Cofu+4Be 


3 


Die Randbedingungen (36) und (37) verlangen, daß die Konstanten A,, B,, C, verschwinden 
und für die Randellipse «, folgende Gleichungen erfüllt sind: 


(69), 
+ &in2u B-2C 
Sie enthalten drei Bedingungen für die Konstanten A, B, C. Die Ausrechnung ergibt 
_ alt —a)e*Coju, 
B= ; 
(4 —a) e" Coju, +8— a) 
2(8—.a) 


Die Spannungsverteilung am Lochrand kann in einfacher Weise angegeben werden; dort ver- 
schwindet o,, mithin kann 0,=0,+0,=2 (1 _ =) AF, oder 


(2) 
= = e"o &of u) co? v+|Sin’u,+ = e % Sin u sin? 
gesetzt werden®).. Nach Umformung ergibt sich 
u 
1 4 — + cos 2» 23). 


p Sin®u, + cos®’v 


2 


9) Bei Unterdrückung des Gliedes mit 43, 


XUM 
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Omax 


x 
03 = . 
Bild 6. Spannungsverlauf in einer gebogenen 
Platte mit Kreisloch. 

Bild 7 (rechts). Spannungsverlauf in einer ; 
gebogenen Platte mit Langloch. 4:9 
Der Spannungsverlauf am Lochrand ist in den mL mh 


Bildern 6 und 7 ersichtlich. Es handelt sich 


um die Fälle (kreisförmiges Loch) und 
Ctgu,=2. Es zeigen sich ähnliche Vertei- 8 
lungen wie beim gelochten Stab. Am Ende 
der kleinen Achse tritt unabhängig vom ? 
Krümmungsmaß stets die gleiche Spannung FG 
1 
1 m 1 05 
02. 


auf, während am Endpunkt der großen Achse « G 
sich — wie zu erwarten war — die Höchst- 


spannung ausbildet, welche mit zunehmender 3 u 


Krümmung zunimmt. Nach Einführung des 


Krümmungshalbmessers o ergibt sich für die 2 
Höchstspannung die einfache Formel 
1 
+.) 7 2 3 4 5 6 7 8 
p 3+ — e Bild 8. Die Höchstspannung in gebogenen Platten mit 
Langloch in Abhängigkeit von der Kerbkrümmung. 
Wieder zeigt sich gegenüber dem entsprechenden Scheibenproblem (gelochter Zugstab ") ) 
eine gemäß dem Faktor Bere geringere Spannungserhöhung. Demnach liegen die Kerb- 
3+— 
m 


spannungen in der gebogenen Platte offenbar stets um 60° (für 1/m =0,) niedriger als 
beim Zugstab. 


VII. Zusammenfassung. Es wird zunächst eine Theorie der dicken Platte aus dem Drei- 
funktionenansatz abgeleitet, welche den Vorteil besitzt, daß die Glieder höherer Ordnung den 
Kraftfluß nicht beeinflussen. Ferner werden für die Randbedingungen, sowie für die Kraft- 
flußintegrale vereinfachte Darstellungen gegeben, mit deren Hilfe die Lösung von Platten- 
problemen gefunden werden kann, ohne auf die in der Regel verwickelten Ausdrücke für die 
Spannungskomponenten eingehen zu müssen. Nach dem so gewonnenen Lösungsverfahren 
werden die beiden Grundprobleme der gekerbten Platte: Gebogene Platte mit beider- 
seitiger Außenkerbe und gebogene Platte mit Langloch einer strengen Lösung 
zugeführt. Die Ergebnisse weisen darauf hin, daß die Spannungserhöhung infolge Kerbwirkung 
bei der gebogenen Platte um etwa 60°) geringer ist als beim Zugstab. Die gewonnenen Form- 
zahlen sind in Diagrammen über dem Krümmungsmaß aufgetragen und bilden die Grundlage 
für die Berechnung der Formzahl für beliebige Abmessungen nach meinem in „Kerbspannungs- 
lehre* gegebenen Interpolationsverfahren. 131 


10) C. E, Inglis; Trans, Instn, Naval Archit., London Bd. 60 (1913), S. 219; Kerbspannungslehre, S. 50. 
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Die experimentelle Ermittlung räumlicher Spannungs- 
zustände an durchsichtigen Modellen mit Hilfe des 
Tyndalleffektes*). 


Von Hermann Josef Menges VDI in Darmstadt. 


A. Theoretischer Teil. 


1. Die Beziehungen zwischen dem Dehnungszustand und dem Indexellipsoid. Bekanntlich 
nimmt ein isotroper, durchsichtiger Körper, wenn er verformt wird, dem Licht gegenüber die 
Eigenschaften eines Kristalles an. Zwischen den Verzerrungsgrößen und dem Indexellipsoid 
bestehen bemerkenswerte Beziehungen, die heute in weitem Maße zur experimentellen Spannungs- 
ermittlung benutzt werden. 


Sind Ay, die drei senkrecht aufeinanderstehenden Hauptbrechungsindizes, so 
lautet die Gleichung des Indexellipsoides, das im Kristall festliegt und auf die Hauptbrechungs- 
achsen bezogen ist: 


Tritt eine ebene Welle linear polarisierten Lichts, deren Wellennormale f sei, in den 
Kristall ein, so wird sie in zwei ebene, linear polarisierte Wellen gespalten, deren Schwingungs- 
ebenen senkrecht aufeinanderstehen und die mit verschiedenen Normalengeschwindigkeiten den 
Kristall durcheilen. Die Abweichungen in den Strahlrichtungen der beiden Wellen dürfen 
nach Untersuchungen Neumannst) für die vorliegende Arbeit außer acht gelassen werden. 
Die Ebene senkrecht zu f durch den Mittelpunkt des Ellipsoids (1) schneidet eine Ellipse aus 
mit den Halbachsen »’ und n’”’ — in Zukunft kurz .„Nebenachsen“ genannt —, die mit den 
Richtungen der dielektrischen Verschiebungen D’ und D” der beiden Wellen übereinstimmen. 
Deshalb sind die aus f und n’ bzw. n’” gebildeten Ebenen die Schwingungsebenen. Außerdem 
bestimmen die Längen der Halbachsen »’ und n” die zugehörigen Normalengeschwindigkeiten 
ce’ und ec”, da allgemein für den Brechungsindex gilt, wenn c, die Lichtgeschwindigkeit im 
Vakuum ist: 


n—=-. 


Es bestehen nun folgende Beziehungen zwischen den Dehnungen und den durch sie 
hervorgerufenen optischen Erscheinungen, die ihren Niederschlag im Indexellipsoid finden: 


1. Die Hauptachsen der Dehnungen &;,, &p, &m Stimmen mit den Richtungen der Haupt- 
achsen des Indexellipsoides überein. 


2. Von allen Fasern, die in der Ebene der obenerwähnten Schnittellipse liegen und vom 
Mittelpunkt ausgehen, weisen zwei aufeinander senkrecht stehende die größte bzw. 
die kleinste Dehnung #’ bzw. &’’ auf, die in Zukunft „Nebendehnungen‘“ genannt seien. 
Die Richtungen von e’ und e’”’ stimmen mit den Richtungen der Nebenachsen überein. 


3. Nach der Neumannschen Theorie gilt die Beziehung: 


K = Stoffkonstante. 


Meßtechnisch besteht demnach die Aufgabe, aus optischen Beobachtungen an einer ebenen 
Welle für einen Punkt im Innern die Schwingungsrichtungen, also die Richtungen der 
Nebenachsen sowie die Differenz n’—n' zu finden. Aus der Durchführung für mehrere 
zweckmäßig gewählte Lichtrichtungen durch einen Punkt sind die Richtungen der Haupt:- 
dehnungen und die Hauptdehnungsdifferenzen zu ermitteln, womit auch die Halbmesser der 
drei Mohrschen Kreise gefunden sind, sofern das Hookesche Gesetz gilt. Es ist demnach 
ausreichend, wenn im folgenden das Verfahren nur für eine Lichtrichtung dargelegt wird. 


*) Nach dem bereits vor längerer Zeit erfolgten Abschluß der theoretischen Betrachtungen zu dieser Arbeit 
wurde mir ein „Letter to the Editor“ von R. Weller, Journal Appl. Physies, Aprilheft Bd. 10 (1939), bekannt, der, 
ohne Einzelheiten zu bringen, auf die Möglichkeit hinweist, den Tyndalleffekt zu Spannungsuntersuchungen zu benutzen. 


ı) F. Neumann: Gesammelte Werke, Bd. 3, Leipzig 1912, 


AUM 


XUM 
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2. Der Polarisationszustand des Lichts und der Tyndalleffekt. 
a) der Polarisationszustand. 


Eine kreiszylindrische 
Scheibe von der Dicke 
Sei einem ebenen 
Spannungszustand un- 
terworfen. Eine engbe- 
grenzte, ebene Welle 
linear polarisierten ein- 
farbigen Lichts durch- 
dringe die Scheibe in 
der Achsenrichtung,also 
normal zur Ebene der 
angreifenden Kräfte 
(Bild 1). Bei der künst- 

lichen Doppelbrechung Fa 
dürfen die dielektri- 
schen Verschiebungen 
und die elektrischen 
Vektoren €’ und €” als 
gleichgerichtet betrach- 
tet werden. Sie sind 

durch die Richtungen ER 

n’ und n” bestimmt. Bild ı. 


Nach dem bekannten Ansatz für ebene fortschreitende Wellen lassen sich die Augen- 
blickswerte e’ und e” der beiden Wellen im Punkte x zur Zeit ? schreiben: 


10) 
cos(ot 
Hierin sind 


2: 
T c„— Kreisfrequenz, Wellenlänge, ce’ = 


—= Normalengeschwindigkeiten. 


‚Ist £ der Winkel, den die Schwingungsebene des einfallenden Lichts mit »’” bildet, und setzt 


man den Absolutwert des elektrischen Vektors gleich Eins, dann sind die Absolutwerte der 
entsprechenden Vektoren der beiden Wellen: E’=sin&, E’=cosE. 
Eine einfache Umrechnung liefert: 


[23 ’ 


Da r im wesentlichen die Zeit enthält und für beide Wellen gleich ist, besteht zwischen 
ihnen die Phasenverschiebung 9. Ihre Verknüpfung mit der Dehnungsdifferenz je Längen- 
einheit ergibt sich unter Beachtung von (2) zu 


d 
K, = Stoffkonstante. 

Aus (3) folgt, daß im Innern eines verformten Körpers das Licht im allgemeinen elliptisch 
polarisiert ist. Der aus beiden Wellen resultierende elektrische Vektor beschreibt während 
einer Periode die Ellipse (3) (in Zukunft kurz „Lichtellipse‘“ genannt), welche auf die Neben- 
achsen bezogen ist. 

Ihre Halbachsen haben bekanntlich die Größen: 


und sind gegen n” um den Winkel y, bezw. y,+ 5 geneigt, der aus der Beziehung folgt 


Jede Lichtellipse ist bestimmt durch die Angabe der Nebenachsen und der Größen und p. 
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b) Der Tyndalleffekt. Fallen z. B. die Sonnenstrahlen auf die feinen Teilchen des 
Rauches einer Zigarre, die an und für sich unsichtbar sind, so werden sie durch Beugung sichtbar. 
Diese Erscheinung nennt man den Tyndalleffekt, dessen Intensität vom Polarisationszustand 
des Lichts abhängt. Angenommen, es sei nur die eine linear polarisierte Welle e’”’ in der 
Scheibe des Bildes 1 vorhanden, die einen Tyndalleffekt zeigen soll. Blickt man senkrecht 
zur Normalenrichtung f unter den Winkel y’ zur Schwingungsebene n”, so weist das in der 
Blickrichtung abgebeugte Licht eine bestimmte Lichtstärke S auf, die verhältnisgleich dem 
Quadrat des elektrischen Vektors E, ist, der dem abgebeugten Licht zukommt. Unter bestimmten 
Voraussetzungen über die Form und Art der beugenden Teilchen läßt sich nachweisen, daß 
zwischen E, und dem Vektor des einfallenden Lichts E’” die Beziehung besteht: 


K = Stoffkonstante. 


Da demnach für E, nur die Komponenten von E” bzw. E’ in Frage kommen, die senkrecht 
zur Blickrichtung stehen, findet man für das Zusammenwirken der beiden Wellen e’ und e”, 
also elliptisch polarisiertes Licht, nach einer einfachen Zwischenrechnung für die Stärke S 
des Tyndalleffektes in der Blickrichtung: 


u=(cos®’2E+ sin?’2 tg2y, —=tg2£&cosp 
Da in Zukunft nur Energieverhältnisse benutzt werden, wurde die auftretende Konstante 


gleich Eins gesetzt. Wird y’=y, bzw. y, + 3 ‚ so nimmt S seinen kleinsten bzw. größten 
Wert an: 


Da aber unter Beachtung von (5a) und (5b) y, mit y, und (7a) mit der kleinen bzw. 
großen Halbachse der Lichtellipse übereinstimmen, erhält man das Ergebnis: 


„Die Blickriehtungen, in denen die kleinste bzw. größte Lichtstärke auftritt, stimmen 
mit den Lagen der großen bzw. kleinen Achse der Lichtellipse überein. Die entsprechenden 
Lichtstärken sind den Halbachsen verhältnisgleich.“ 


Durch Messungen der Stärke bzw. der 
Menge des abgebeugten Lichts kann man dem- 
nach grundsätzlich die Lichtellipse sowohl ihrer 
Größe als auch ihrer Lage nach in jedem 
Punkt des Körpers ermitteln. 

Es ist zweckmäßig, den experimentell be- 
stimmbaren Wert u als Betrag einer im Körper 
festliegenden gerichteten Größe u zu betrachten, 
die in Richtung der großen Achse der Licht- 
ellipse weist. Gegenüber einer einmal fest- 
gelegten, aber sonst beliebigen Bezugsebene 
bzw. Grundrichtung im Körper besitzt u den 
Winkel &, =ö-+y,, wenn ö die Neigung der 
Achse n’’ gegen die Grundrichtung ist. 

Durch Verdoppelung aller Winkel geht 
u in u* über, ohne seinen Betrag zu ändern. 
Unter Beachtung der in (7) ausgedrückten 
Beziehung für « erhält man Bild 2, das für 
die spätere Auswertung grundlegend ist. Bild 2. 


3. Der räumliche Spannungszustand. Es sei angenommen, daß der zylindrische Körper 
- nach Bild 1 nun einem allgemeinen räumlichen Spannungszustand unterworfen sei. Dann 
ändern die Dehnungsellipsoide im Körper ihre Lage und Größe von Punkt zu Punkt. Folglich 
sind auch längs des Lichtstrahles die Nebenachsen zweier um Ax auseinanderliegender 
Nachbarelemente x, und &, um den Winkel Aö gedreht. Ferner erhalten die beiden Licht- 
wellen beim Durcheilen des Elements Ax eine zusätzliche Phasenverschiebung A, die nur 
von der Differenz der Nebendehnungen €’ — e” abhängig ist und sich von Ort zu Ort ändert. 


XUM 
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Die Aufgabe besteht darin, die Funktionen 
d 


welche den Spannungszustand längs des Lichtstrahles bestimmen, experimentell zu ermitteln. 
Da durch Messung nur u als Funktion von x bzw. ö bekannt ist: 
u=f,(&) (9) 
muß man versuchen, (8) aus (9) zu gewinnen. Man beachte, daß jetzt die gesamte Phasen- 
verschiebung, die in Zukunft nicht mehr @, sondern 2 genannt sein soll, sich aus 2 Anteilen 
zusammensetzt: nämlich einem Teil, der auf die Drehung der Nebenachsen zurückzuführen 
ist, und einem zweiten, der aus den Differenzen der Nebendehnungen folgt. In Bild 3 sind 
die Verhältnisse an zwei benachbarten Elementarscheiben 1 und 2 von der Dicke Ax dar- 
gestellt. Die Lichtellipse am Ende der ersten Scheibe ist bestimmt durch die Größen £&, und 
2, bzw. u,. Da das Licht beim Eintritt in die zweite Scheibe nur in Richtung der Neben- 
achsen schwingen kann, muß man die Ellipse, die ihre Gestalt und ihre Lage im Körper 
ebenso wie u, beibehält, auf die neuen Achsen beziehen, d. h. es ändern sich £, in &, und 
2, in Q,. Beim nun folgenden Durcheilen der Scheibe 2 ändert sich £, nicht, da die Neben- 
achsen konstant sind. Infolge der auftretenden Differenz der Nebendehnungen besteht jedoch 
am Ende dieser Scheibe eine neue Phasenverschiebung 2,= Q2,'+ Ap,so daß das Licht nach 
einer neuen Ellipse mit einem neuen u, schwingt, die durch £, und 2, (in Bild 3 zufällig 
a/2), bestimmt ist. 


Ende der 1. Scheibe Anfang der 2.Scheibe Ende der 2.Scheibe 


In Bild 4 ist das in Bild 2 bereits 
dargestellte Diagramm auf diese Verhält- 
nisse erweitert worden. Das rechtwinkelige 
Dreieck PC, D, mit &, und 2, gilt für das 
Ende der ersten Scheibe, ist also auf die 
Nebenachsen n,' n,” bezogen. Um die 
Bestimmungsgrößen &, und 2, von u, in 
bezug auf die Nebenachsen n,’' n,", also 
am Anfang der zweiten Scheibe zu gewinnen, 
projiziere man u,* auf die unter 2, zur 
Grundrichtung geneigte Grade (n,”), dann 
ist PB,=2E£,, C,B,=sin 2£, und 
0‚„D,=sin2£,cos 2,, da die in (7) aus- 
gedrückte und in Bild 2 dargestellte Be- 
ziehung für jedes Nebenachsensystem gelten 
muß. Diese Größen sind dem Diagramm 
bei gegebenen u,* und (n,”) zu entnehmen, 

Da nun £, für die zweite Scheibe 
konstant bleibt, kann der Vektor u,*, der zum 
Ende der zweiten Scheibe gehört, mit seiner 
Spitze D, nur auf der Graden C, B, liegen. Bild 4. 


& 
> 
n," Nur 
Bild 3. 
\ 
® 
— | 
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Damit ist aber das für die Bestimmung der Nebenachsen sehr wichtige Ergebnis gewonnen: 
„Das Lot von P auf die Differenz Au* zweier benachbarter Vektoren u,* und 
u,*:du*=u,* — u,* bildet mit der Grundrichtung einen Winkel, dessen Halbierende die zu 
u, gehörige Nebenachse n," ist.“ 
302 Ist u eine stetige Funktion von ö, so bilden 
die Spitzen der Vektoren u* eine Kurve, kurz 
x : „Spitzenkurve‘“ genannt, und 4 u* geht in du* über. 
Dann lautet der oben ausgesprochene Satz: 

„Das Lot von P auf die Tangente an die 
Spitzenkurve im Punkt u* bildet mit der Grund- 
richtung einen Winkel, dessen Halbierende die zu 
u gehörige Nebenachse n’ ist.“ 

Es seien noch zwei Differentialgleichungen 
abgeleitet, die für das optische Geschehen in einem 
unter räumlicher Spannung stehenden Körper grund- 
legend sind. 

In Bild 5 sind im wesentlichen wieder die 
Diagramme zweier benachbarter Vektoren u,* und 
u,* und die Spitzenkurve F dargestellt, die zu zwei um 
46 cos24ö-= 1) gegeneinander 
geneigten Achsensystemen gehören. Man kann 


‚=UG,E+ED, 


Bild 5. ausdrücken durch: 


sin28,c08 Q,+Aysin . . . . (10), 
da 2’ = Q,— Ag ist. 


Man ersetze &, durch £&, unter Vernachlässigung der Produkte kleiner Größen 


cos2&, (,E=cos2£, — sin2£,cos Q,-2 46 I. 
sin2&,=sin2&,+246-c0s2£, cos 2, 
Eine einfache Zwischenrechnung liefert, wenn man (11) in (10) einführt: 
Da nun = 2, +42 ist, folgt: 
13). 
- dö 


Gl. (12) zeigt deutlich, daß der Zuwachs an Phasenverschiebung d 2, den das Licht beim 

Durcheilen eines Elementes erhält, nur zu einem Teil dp durch die Differenz der Neben- 

dehnungen d. h. die Größe des Spannungsellipsoides bedingt ist. Der mit dö behaftete 

Ausdruck stellt den Zuwachs dar, der aus der Drehung der Spannungsellipsoide folgt. Aus 

diesem Grund lassen sich die Verfahren der ebenen Spannungsoptik (dö—=0), z. B. die 

Kompensation der Phasenverschiebung, nicht ohne weiteres auf räumliche Zustände übertragen. 
Beachtet man, daß ,=&, +48 ist, so geht die zweite der Gl. (11) über in: 


dE=— cos Q,- 
oder — cos dd. 
Die Integration liefert mit der Konstanten n: 


0 
Die beiden Gl. (11) und (13), die als die Grundgleichungen der Optik in einem 
verspannten Körper betrachtet werden dürfen, sind bereits von Neumann’) auf anderem 
Weg und in etwas anderer Form entwickelt worden. 


4. Die Auswertung. Hat man aus Messungen die Vektoren u längs des Lichtstrahls 
in einer ausreichenden Zahl von Punkten x gefunden, so zeichnet man unter Beachtung der 
Grundrichtung im Modell das Polardiagramm der u*, wovon ein Ausschnitt etwa das Bild 4 


?2) F. Neumann: Gesammelte Werke, Bd. 3, Leipzig 1912. 
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bringt. Ferner trage man den Einheitskreis ein. Zu zwei benachbarten Punkten x und + 4x 
des Körpers gehören die beiden Vektoren u,* und u,*. Man verbinde ihre Spitzen D, D, und 
fälle die Senkrechte PC,. Die Anwendung des Satzes des Abschnpittes A, führt zu dem 
zugehörigen mittleren Kreuz der Nebenachsen. 

Damit ist der erste Teil der in A, gestellten Aufgabe (8) erfüllt. Der zweite Teil, 


nämlich zr- f.(x) zu finden, läßt sich mit Hilfe des Diagramms von Bild 4 ebenfalls rasch 
lösen. Es bestehen die Beziehungen: 


C,„B,=sin2£,, cosQ'= 
C, D,=sin 2£, cos Q,, (AB). 
C,B, 


Die Strecken C,B,, €, D, und C,D, entnimmt man dem Diagramm und rechnet nach 

(14) und (15) 2,’ und 2, aus, wobei man nur noch auf das Vorzeichen zu achten hat. Da 

nun Ago=Q,— (2, ist, hat man damit auch A gefunden. Das zugehörige Intervall Ax ist 
RK, 


ebenfalls bekannt, also liefert Formel (5): 2? es (€ —e”) für jeden Punkt die Differenz 


der Dehnungen. 

Die Lage der Schwingungsebene, oder allgemeiner der Polarisationszustand des Lichtes, 
das in den Körper eindringt, tritt in diesen Darlegungen nirgends explizit auf. Man kann 
jedoch nachweisen, daß die Integrationskonstante n der Gl. (13) im wesentlichen die 
Lage der Schwingungsebene zur Grundrichtung im Körper enthält. In u ist demnach die 
Stellung des Polarisators enthalten. Wie im folgenden bei der Besprechung der experimen- 
tellen Ergebnisse auch gezeigt wird, liefern verschiedene Schwingungsebenen des einfallenden 
Lichts auch verschiedene Funktionen u—=f(,) und damit andere Spitzenkurven. Da aber 
die Nebenachsen und die Differenzen der Nebendehnungen physikalisch gegeben sind, müssen die 
Achsen (ö) und die Größe Ag invariant in bezug auf den Polarisationszustand des einfallenden 
Lichtes sein. Daraus folgen eine Reihe von Gesetzmäßigkeiten über die Beziehungen zwischen 
u und diesem Zustand, von denen man Gebrauch machen kann. So müssen z. B. die Tangenten an 
entsprechenden Punkten verschiedener Spitzenkurven, die mit verschiedenem polarisiertem Licht 
aufgenommen wurden, parallel sein. Dies ist eine gute und auch schon äußerlich leicht fest- 
zustellende Kontrolle für die Richtigkeit einer Messung (s. Bild 8). 


B. Experimenteller Teil. 

1. Beschreibung der Meßanordnung: Die Versuche wurden mit einer Anordnung gemacht, 
die aus vorhandenen Teilen zusammengestellt werden mußte. Aus diesem Grund sind die 
Möglichkeiten, was die Genauigkeit und Schnelligkeit der Messungen anbelangt, bei weitem 
nicht ausgeschöpft. Die Messungen sollen lediglich den Beweis erbringen, daß grundsätzlich 
keine meßtechnischen Schwierigkeiten die Durchführbarkeit des Verfahrens unmöglich machen. 

Als Lichtquelle @ (Bild 6) wurde 


eine Quecksilberdampflampe (Nieder- £ 

druck) älterer Ausführung benutzt. Die + 
Verwendung einer solchen Lampe, deren I 
Licht bekanntlich starke kurzwellige, Ah 
blaue und violette Linien aufweist, gibt ; 


deshalb eine gute Ausbeute, da die Kon- 
stante der Formel (6), d. h. aber die ab- 
gebeugte Lichtmenge dem Kehrwert der 


vierten Potenz der Wellenlänge ver- N 
Das Licht ging durch ein Filter ® | 
(F), das im wesentlichen nur Licht der ! 1 
Wellenlänge A = 436 m u durchließ. 

Ein Fotoobjektiv A mit verhältnis- 
mäßig geringer Lichtstärke, aber großer Bild 6. 
Brennweite, sammelte das Licht und ent- 
warf auf der Kreisblende B von 2 mm Durchmesser ein Bild des Lichtbogens der Lampe. 
Unmittelbar vor der Blende befand sich ein Nicolsches Prisma, dessen Schwingungsebene 
der Blende gegenüber verdreht werden konnte. Das polarisierte Licht ging durch das un- 
mittelbar hinter der Blende befindliche Modell M, das einem Spannungszustand unterworfen war. 
Nicol-Blende und Modell konnten um den Lichtstrahl als Achse gemeinsam gedreht werden. 
Der Winkel $ der Grundrichtung gegen eine feste Marke war an einer Gradteilung ablesbar. 
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Als Modelle fanden gegossene Trolonzylinder von 3 cın Durchmesser und 3 cm Länge 
Verwendung. Um für diese ersten Versuche den Bau einer Belastungsvorrichtung zu um- 
gehen, wurde der Eigenspannungszustand geprüft, mit dem jede Probe bei der Anlieferung 
behaftet ist. Es ist kaum ein Nachteil, daß dieser Zustand keine unmittelbare Nachprüfung 
gestattet, da die Richtigkeit der theoretischen Untersuchungen nach der unter A4 dar- 
gelegten Weise bei verschiedenen Lagen der Schwingungsebenen des Polarisators nach- 
zuprüfen ist. 


Im rechten Winkel zum Lichtstrahl befand sich das Mikroskop Mi, das unter Zwischen- 
schaltung einer geeigneten Blende ein kurzes Stückchen des Lichtstrahles auf der lichtempfind- 
lichen Schicht der Photozelle Ph abbildete. Mikroskop und Photozelle waren längs des 
Lichtstrahles beweglich. Als Photozelle fand eine hochisolierende Kalium-Vakuum-Zelle mit 
ausreichender Empfindlichkeit im Blauen Verwendung. Leider war es dem Verfass»r nicht 
möglich, geeignete Hochohmwiderstände zu erhalten. Deshalb wurde die abgebeugte Licht- 
menge nach dem Aufladeverfahren mit Hilfe eines Lindemann-Elektrometers E gemessen. 
Der positive Pol der Saugbatterie war mit der Anode der Photozelle verbunden, der negative 
Pol und das Elektrometergehäuse lagen an Erde. Die Kathode führte zur Nadel des Elektro- 
meters. Unterbricht man ihre Erdung a, so fließt ein Aufladestrom über die Zelle, der ein 
Anwachsen der Potentialdifferenz zwischen Nadel und Erde bewirkt. Die Zeit, welche der 
Zeiger des Elektrometers benötigt, um eine bestimmte Anzahl Skalenteile zu durchlaufen, ist 
verhältnisgleich dem Widerstand der Photozelle und damit umgekehrt verhältnisgleich der 
abgebeugten Lichtmenge. Alle Teile der Anordnung waren, soweit erforderlich, lichtdicht ein- 
geschlossen und elektrisch abgeschirmt. 


2. Durchführung und Ergebnis der Messungen. Den bisherigen Betrachtungen lag ein 
„idealer“ Tyndalleffekt zugrunde: 


1 
—2y,)). 


"Für linear polarisiertes Licht 2= 0 müßte S für eine bestimmte Lage des Modells zur Grund- 
richtung verschwinden. In Wirklichkeit ist dies nicht der Fall, da infolge von Streulicht 
oder nicht völliger Polarisation usw. noch eine kleine Lichtmenge S, — oder genauer gesagt — 
noch ein entsprechender kleiner Anzeigewert vorhanden ist. Es hat sich ergeben, daß die 
Verhältnisse gut erfaßt werden durch die Formel: 


S=8,+48,(1- 


u= (cos®’28+sin?2 


Außerdem war es mit den zur Verfügung stehenden Mitteln nicht möglich, die Lichtmenge 
der Quecksilberlampe längere Zeit völlig konstant zu halten. Es wurde folgender Weg be- 
schritten: 


Eingehende Vorversuche zeigten, daß bei einer Änderung der Helligkeit der Lampe 
oder der Empfindlichkeit des Elektrometers, sich sowohl S als auch S, und S, mit denselben 
Verhältnisziffern änderten. Man kann also einer geeignet gewählten Lichtmenge willkürlich 
den Wert 1 beilegen und alle anderen Messungen auf sie umrechnen. Um eine solche Ver- 
gleichsgrundlage zu erhalten, wurde eine mit dem zu prüfenden Modell völlig gleichartige 
Trolonprobe etwa 2 Stunden auf 60° bis 70° Wärme gehalten und dann sehr langsam abge- 
kühlt. Eine Nachprüfung ergab, daß sie praktisch frei von Eigenspannungen war. Für diesen 


Be ı Zylinder wurde an verschiedenen Stellen 
unverformtes x die Funktion (16) S=f(y') bzw. 
S=f(ß), die sich übrigens praktisch 
as deckten, durch Drehen des Modells für 
eine ausreichende Zahl Winkelstellungen 
46 ip experimentell bestimmt (s. Bild 7). 
05: Per Maßstab der Ordinate S wurde so 
“IS r (Pınkt3Kurved) gewählt, daß die doppelte Amplitude der 
3  Cosinusfunktion, d.h. u für linear polari- 
a | siertes Licht gleich Eins gesetzt wurde. 
vo Als Umrechnungsziffer für die Messungen 


am verformten Modell diente das Ver-. 
Bild 7. hältnis der ÖOrdinaten der Mittellinie 
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1 
(+3 5,) zur entsprechenden Ordinate gemessen am unverspannten Trolonzylinder. Um die 


Meßergebnisse auszugleichen, wurde jeweils die beobachtete Kurve durch Ausplanimetrieren 
in die flächengleiche trigonometrische Funktion umgewandelt. Es zeigten sich dabei nur 
geringfügige Abweichungen von dem nach (16) geforderten Verlauf. Wie aus Bild 7 ersichtlich, 
ist auch der Winkel /,, der die Richtung des Vektors « festlegt, leicht zu ermitteln. 


Bild 8 enthält die Ergebnisse der Messungen für eine Meßstrecke #,—4mm bei drei 
verschiedenen Nicolstellungen in Gestalt der Funktionen «*—f(ß,). Man erkennt rein 


Zahlentafel. 


Zr; Meß- | Meß- | Meß- | Meß- 
reihe | reihe [reihe reihe 


6° mm 


85,25 87,7° | 13,8 | 14,5 
77,25 | 795 | ı29 133 


1mm 
78 78 11,5 | 11,2 


785 | 75 | ı28| 139 


Bild 8 (links). 


äußerlich, daß die Tangenten in entsprechenden Punkten der 3 Kurven gut übereinstimmen. 
Dabei sei darauf hingewiesen, daß die Meßreihe © mehrere Wochen vor den Reihen A und B 
aufgenommen wurde. Dadurch erklären sich kleinere Abweichungen, da infolge von Aus- 
gleichserscheinungen die Spannungszustände sich nicht mehr voll decken. In der Zahlen- 
tafel 1 sind die aus den beiden Kurven A und B berechneten mittleren Werte für die 
Winkel ö (Lage der Nebenachsen des Spannungszustandes) und für die Größen er angegeben. 
Mit Rücksicht auf den nur für einen ersten Versuch bestimmten Aufbau der Anordnung sind 
die Ergebnisse als durchaus befriedigend und für die Zukunft recht erfolgversprechend zu 
bezeichnen. 


C. Zusammenfassung. 


Die vorliegende Abhandlung befaßt sich mit der Anwendung des Tyndalleffektes zur 
Ermittelung räumlicher Spannungszustände Nach einer kurzen Einführung in die grund- 
legenden Gesetze wird ein einfaches, zeichnerisches Verfahren entwickelt, um die Meßergeb- 
nisse auszuwerten. Anschließend folgt eine Beschreibung der Versuchsanordnung und die 
Besprechung der ersten so gewonnenen Ergebnisse. 

Wie schon betont, lassen sich die Genauigkeiten und die Schnelligkeit der Messungen 
ganz wesentlich erhöhen, da man bei der Benutzung neuzeitlicher Lichtquellen und einer 
besseren Optik eine wesentliche Erhöhung der abgebeugten Lichtmengen errechnen kann. 
Durch Verwendung von Hochohmwiderständen in Verbindung mit dem Elektrometer oder auch 
von geeigneten Verstärkern und Photozellen mit höherer Empfindlichkeit läßt sich die Zeit 
zur Durchführung einer Messung durch Anwendung einer automatischen Aufzeichnung auf 
wenige Minuten herabdrücken. 

Es ist möglich, auch andere Auswerteverfahren zu benutzen. So kann man, ähnlich wie 
in der ebenen Spannungsoptik, die auch beim Tyndalleffekt auftretenden Interferenzfarben 
bei Anwendung von weißem Licht heranziehen. 

Geeignet erscheint das Verfahren auch in Verbindung mit dem von Oppel?) angegebenen 
„Erstarrungsverfahren“ zu sein, da sich zum Auswerten ein Zerschneiden des Modells erübrigt. 

Bei beliebig geformten Modellen muß man, um Brechungen zu vermeiden, wie bekannt, 
die Proben in Zylinder einschließen, welche mit einer Flüssigkeit gefüllt sind, die gleichen 
Brechungsindex wie der Modellstoff hat. 186 


3) Oppel, Forsch. Ing.-Wes. Bd. 7 (1936), S. 270. 
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Von den freien Schwingungen eines Kreiselpendels 
bei endlichen Ausschlägen*). 


Von U. T. Bödewadt in Göttingen. 


$ 1. Einleitung. Als Kreiselpendel bezeichnet man einen (symmetrischen) Kreisel, dessen 
Achse durch quasielastische Kräfte (Schwerkraftkomponenten) an eine Gleichgewichtslage ge- 
bunden ist, um welche sie Schwingungen ausführen kann. Dieser Fall kommt in der Technik 
mehrfach vor, z. B. in Gestalt des von Schlick erfundenen Stabilisierungskreisels, der die 
Drehschwingungen eines Schiffes um seine Längsachse, das Rollen oder Schlingern, dämpfen 
soll. Die Bewegung der Achse eines solchen Schiffskreisels kann man durch Angabe von 
zwei Winkeln y, ® in zwei zueinander senkrechten Ebenen darstellen, deren Schnittgerade 
y=®d=0 die Ruhelage der Kreiselachse ist; sie gehorchen den Differentialgleichungen ') 


d? 
+ OH sin y— N 
(1.1). 


.d» 
sind N (ein y)=0 


Darin bedeuten: y den Rollwinkel des Schiffes, 9 den Pendelungswinkel des Kreisels; J, j die 
zugehörigen Trägheitsmomente des Schiffes bzw. Kreisels; Q die Wasserverdrängung (Masse), 
H die metazentrische Höhe des Schiffes, q und h die entsprechenden Größen beim Kreisel; 
N das Drehmoment des Kreisels, d. i. das Produkt aus seiner Umdrehungszahl f und seinem 
Trägheitsmoment T um die Rotationsachse. 


Dieses Gleichungssystem ist bisher nur für kleine Winkel y, # untersucht worden, wenn 
der Sinus durch den Bogen ersetzt werden konnte, wodurch (1.1) überging in 


d’y d 
Jet 


(1.2). 
d 
IHN 


Wenn man nun auch für größere Winkel die genaueren Gl. (1.1) durch das System (1.2) er- 
setzte, so war man sich wohl der in dieser Vereinfachung beschlossenen Vernachlässigung 
bewußt — unterscheiden sich doch z. B. die letzten Glieder, die Kopplungsglieder, in (1.1) von 
denen in (1.2) um den Faktor cos # bzw. cosy, so daß im Durchschnitt der Kopplungsgrad 
gegenüber den Werten für kleine Schwingungsweiten herabgesetzt ist. Wenn nun auch zu 
vermuten ist, daß diese Vernachlässigung nicht sehr ins Gewicht fällt’), da bereits der An- 
satz (1.1) vereinfachende physikalische Annahmen enthält, so ist es doch von Interesse, den 
Einfluß der mathematischen Vernachlässigung festzustellen, um die Zulässigkeit der physika- 
lischen Vereinfachungen gesondert von der der mathematischen Vereinfachungen beurteilen 
zu können. Als ersten Schritt in dieser Richtung habe ich es daher im folgenden unter- 
nommen, das nicht-lineare System (1.1) auf seine Eigenschwingungen zu untersuchen; die 
nächsten Schritte wären Einführung einer Reibungsdämpfung und Behandlung der erzwungenen 
Schwingungen. Die Anregung zu dieser Untersuchung wurde mir von Herrn Prof. Dr. Schuler 
zuteil, dem ich auch für seine fördernde Kritik während der Arbeit zu Dank verpflichtet bin, 
ebenso wie Herrn Prof. Dr. Kaluza für das mir stets bewiesene Interesse und für seine wert- 
vollen Ratschläge und Hinweise. 

Entsprechend den Verhältnissen in der Praxis, wo zwar die Kreiselausschläge # groß 
sind, die Schiffsschwingungen aber noch klein bleiben, beziehe ich mich statt auf (1.1) auf 
das vereinfachte System 


d’ 


d (1.3); 


*) Diese Abhandlung wurde von der Mathematisch-Naturwissenschaftlichen Fakultät der Universität Göttingen 
als Habilitationsschrift angenommen. 

1) F. Klein und A. Sommerfeld: Über die Theorie des Kreisels. Leipzig 1914, S. 800. 

2) A. a. O.1), S. 833 bis 834. 


| 


2. . Math. Mech. > 
Ba. Nr. er Bödewadt, Von den freien Schwingungen eines Kreiselpendels 219 


doch ist die im folgenden angewandte Methode nicht nur ohne weiteres auf das System (1.1) 
übertragbar°), sondern auch auf andere nicht-lineare Systeme gekoppelter Schwingungen. 


$ 2. Der lineare Grenzfall. Wenn man in (1.3) neue Veränderliche und Konstanten ein- 
führt durch 


V3=? 
so gewinnt man die Form?) 
— ksny—0, 
jy+bsiny+k&—0 


Einer Verkürzung des Zeitmaßes entspricht eine Verkleinerung von a, b, k,c im selben Maße, 
Es kommt daher im wesentlichen nur auf das gegenseitige Verhältnis und nicht auf die ab- 
solute Größe dieser Koeffizienten an. Im linearen Grenzfall y>0 wird aus (2.2) das dem 
System (1.2) gleichwertige Differentialgleichungssystem 


— ky=0, 
(2.3), 
dessen Lösungen gegeben sind durch 
(2.4) 
mit beliebigen Konstanten M,, M,, t,, t,. Darin ist 


und die Frequenzen ®,,, ®,, Sind die positiven Wurzeln der Gleichung 


Diese im Grenzfall kleiner Schwingungen (y—0) sich einstellenden Frequenzen &,,, ®,, sollen 
weiterhin die „Grenzfrequenzen“ heißen. — Bei passenden Anfangsbedingungen kann eine der 
beiden Konstanten M,, M, verschwinden; es gibt also zwei verschiedene periodische Lösungen 
für das System (2.3), jede noch mit willkürlicher Phasenkonstante: ?, bzw. t,, und willkürlicher 
Amplitude: M, bzw. M, in y(f), oder M,r, bzw. M,r, in «(f). Die allgemeine Lösung kann 
aus diesen zwei periodischen Lösungen linear zusammengesetzt werden, ist also in der Regel 
(nämlich, wenn das Verhältnis der Grenzfrequenzen irrational ist) eine fastperiodische Funktion 
der Zeit £ — und zwar eine der einfachsten Art, weil in ihrer Fourierreihe nur zwei Fre- 
quenzen vorkommen. Haben die Grenzfrequenzen ein rationales Verhältnis, so ist die allge- 
meine Lösung (2.4) insbesondere reinperiodisch. 

Die langsamere der beiden Schwingungen soll hier wie üblich als Präzession bezeichnet 
werden °), ihre Frequenz mit die schnellere heiße Nutation°), ihre Frequenz sei Nach 
(2.6) ist also 


1 en 


3) S. u. $ 10. 

4) Die Punkte bedeuten wie üblich Ableitungen nach der Zeit t. Die Überstreichung von sin y soll eine 
Klammer ersetzen. 

5) A. a. 0.1), S. 803, werden diese beiden Schwingungen als „Hauptschwingung“ und „Nebenschwingung‘“ unter- 
sehieden. — Es ist bei dieser Festsetzung durchaus möglich, daß bei größeren Ausschlägen etwa die hier als „Präzession* 
bezeichnete periodische Lösung über der Grundsehwingung starke Oberschwingungen von einer Frequenz enthält, die 
der Größenordnung nach mit der Grundfrequenz der „Nutation* übereinstimmt; da aber diese Oberschwingungen nicht 
selbständig vorkommen, sondern nur als Bestandteil der gesamten „Präzessions“-Schwingung, so soll diese gesamte 
Schwingung als reine Präzessionsschwingung ohne Nutationskomponente angesprochen werden, mag auch dem Augen- 
schein nach die Unterscheidung schwerfallen. 
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Daraus lassen sich die Ungleichungen entnehmen: 
0< < Min b?), | 


. (283). 
Max (a?, b°)< 


Die beiden demnach stets reellen Kopplungsfrequenzen &,,, ®,, Schließen ihre Grenzwerte 
für k>0, die „ungekoppelten“ Frequenzen a und b, zwischen sich ein; mit wachsendem k 
rücken sie immer weiter auseinander. 

In dem bis jetzt betrachteten Grenzfall y>0 läßt sich alles Wesentliche über die freien 
Schwingungen des Systems (2.3) bereits aussagen, wenn man die beiden periodischen Lö- 
sungen (2.5) kennt, da diese sich in jeder Lösung (2.4) einfach überlagern. Da im System (2.2) 
die Funktion sin y nicht linear in y ist, so kann hier eine lineare Zusammensetzung der all- 
gemeinen Lösung aus den etwaigen periodischen Lösungen nicht mehr in Frage kommen. 
Dennoch wird auch für das nichtlineare System (2.2) die Kenntnis periodischer Lösungen be- 
reits vieles zu überblicken gestatten. 


$ 3. Der Ansatz für die periodischen Lösungen. Es sollen jetzt diejenigen streng perio- 
dischen Lösungen des Systems (2.2) gesucht werden, die für y>0 in Form und Periode stetig 
in eine der beiden periodischen Lösungen (2.5) des Systems (2.3) übergehen®) — die Existenz 
solcher Lösungen, die physikalisch zu vermuten ist, zunächst vorausgesetzt. 

Da die Gl. (2.2) die Zeit £ nicht explizit enthalten, so ist in der Lösung eine willkürliche 
Phasenkonstante; diese werde so gewählt, daß zur Zeit ?=0 die Kurve sich gerade im Um- 
kehrpunkt befindet. Ferner soll die Lösung dieselbe Symmetrie haben wie die entsprechende 
Lösung des Systems (2.3), d.h. wie die Kosinuskurve’): also periodisch mit der Frequenz o, 
symmetrisch zur Geraden t=0 und, bis auf Verschiebung um eine halbe Periode, auch zur 
Geraden y—=0. Der Ausschlag im Umkehrpunkt soll in Bogenmaß mit M bezeichnet werden °). 


Die Randbedingungen sind also 


Nach (2.4), (2.5) folgt aus „>0, d.h. M>0, auch <>0; x und y werden also den Faktor M 
enthalten. Daher führe man, um für M>0 eine endliche Vergleichsfunktion zu haben, neue 
Veränderliche ein durch 
z=Mu; 
Man setze ferner 
siny sinMz 
= (M >20), 
M 
S(2)=2 (M=0) 


Dann erhält man aus (2.2) das System 


S(e)+kü—0 84), 


M>0 


6) Durch diese Forderung des stetigen Zusammenhanges mit den kleinen Schwingungen um die Ruhelage y = 0 
beschränkt man sich auf die — technisch viel wiehtigeren — sich nieht überschlagenden Schwingungen mit |y|<z, 
da diese mit den sich überschlagenden Schwingungen wohl nicht stetig zusammenhängen, sondern in gewissem Sinne 
eine gesonderte Lösungsmannigfaltigkeit bilden. Dies ist jedenfalls für die erzwungenen Sehwingungen bei fehlender 
Kopplung der Fall, vgl. darüber R. Iglisch: Die erste Resonanzkurve beim Duffingschen Schwingungsproblem. 
Math. Ann. Bd. 112 (1956), S. 221 bis 246. 

?) Die Fourierreihe einer solchen periodischen Lösung kann nur die Kosinus der ungeraden Oberschwingungen 
zur Grundfrequenz ® enthalten. Es mag scheinen, als giugen durch diese Symmetrieforderung gewisse periodische 
Lösungen, etwa solche der Form ce, coswt+cgcos2wt+... verloren. Indessen gibt es keine periodischen Lösungen 
mit geraden Oberschwingungen, wenn außerdem > wo verlangt wird — s. u. $ 8. 


s) Wo M in Graden angegeben wird, mag das Zeichen ® als Abkürzung für den Faktor en stehen. 


XUM 
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Die Gl. (3.4) lassen sich auch als ein System von Differentialgleichungen 1. Ordnung 
schreiben; man setze etwa u=x,, 2=x, und hat 
2, =f,(2,, 2, 2, 2)=—- 
=f, (2,20, 2, 


(3.6). 


Nun sind nach (3.3) S(z) und damit die rechten Seiten f,,...f, von (3.6) ganze analytische 
Funktionen ihrer Argumente und des Parameters M?; man erhält für S(z) aus (3.3) die Reihen- 
entwicklung 


3 5 7 


Zu jedem endlichen Wert von A gibt es daher durch jeden Punkt (x, (f,), (f,), (to), &, (t,)) 
eine analytische Lösung, die auch vom Parameter A analytisch abhängt. Mit anderen Worten: 
es gibt im Endlichen keinen (auch keinen komplexen) Wert von 4, der für die Lösungen von 
(3.6) und damit von (3.4) eine singuläre Stelle wäre — jede Lösung von (3.4) kann daher in 
eine beständig konvergente Reihe nach wachsenden Potenzen von 4 entwickelt werden. Die 
Grenze A=n°, d.h. M=180°=n, jenseits derer keine durch die Ruhelage y„=0 gehenden 
oszillatorischen Lösungen zu erwarten sind, wird sich — nach den Verhältnissen beim gewöhn- 
lichen Pendel zu urteilen — nicht durch Existenzforderungen, sondern durch Realitätsforde- 
rungen geltend machen. 

Für A=0 läßt das dann mit (2.3) gleichbedeutende System (3.6) die periodischen 
Lösungen (2.5) zu. Nach dem Stetigkeitssatze von Poincare°) existieren daher in der Um- 
gebung von Ä=0 periodische Lösungen mit i. A. von ®, verschiedener Frequenz &® =» (}), 
die für A>0 stetig in die periodischen Lösungen für A=0 übergehen. Diese periodischen 
Lösungen wie auch die Frequenz » lassen sich diesem Satz zufolge in absolut konvergente 
Potenzreihen nach dem Parameter entwickeln, solange er genügend klein bleibt; über den 
Konvergenzbereich wird Näheres nicht ausgesagt '"). 


Aus den Gl. (3.4) kann man die Veränderliche « eliminieren, indem man die erste der 


beiden Gleichungen nach ? differenziert und mittels der zweiten Gleichung «& und & durch z 
ausdrückt; übrig bleibt eine Differentialgleichung 4. Ordnung für z: 


tab ......... 88. 
Darin ist S(z) die durch (3.3) erklärte, durch (3.7) dargestellte, von dem Parameter A abhängige 
Funktion. Die für hinreichend kleine Parameterwerte existierende periodische Lösung mit 
den Randbedingungen (3.5) hat die Form 


Auch für die Frequenz hat man eine solche Reihe anzusetzen: 


der erste Koeffizient dieser Reihe, &,, muß wegen (3.5) mit dem durch (2.6) gegebenen identisch 
sein und wurde deshalb schon ebenso bezeichnet. 


Setzt man die Reihe (3.9) in (3.7) ein und ordnet nach den Potenzen von 4, so erhält man 
. 8.11). 
2% 2 o-ı 24 5040 
Über die Konvergenz dieser Reihen läßt sich ohne nähere Untersuchung nur noch 
folgendes sagen: Die Reihe (3.7) ist für alle endlichen Werte von A konvergent. Wenn die 
Lösung B. 9), (3.10) nicht für alle A< rn? (also alle M< 180°) konvergiert und die Bedingungen 


6 
+ 


») H. Poincars: Les methodes nouvelles de la mecanique eeleste. Paris 1892, Bd. I, S. 89 bis 95. 

10) Die neueren Untersuchungen von Painleve6, Popoff (z. B. Math. Ann. Bd. 101 (1929), S. 570 bis 578) u. a. 
riehten sich auf die Natur singulärer Stellen der Lösung, wenn als Randbedingungen die Werte für.t=0 gegeben sind 
— nicht bei der Randbedingung der Periodizität. 
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(3.5) erfüllt, so gibt es einen kleineren Wert 4*, jenseits dessen keine Konvergenz stattfindet. 
Dann existieren überhaupt keine Lösungen der verlangten Art für 4>4*, auch nicht in anderer 
Form als der einer Potenzreihe nach A. Denn jede solche Lösung hat eine bestimmte end- 


liche Periode en erfüllt also bestimmte endliche Randbedingungen (3.5) und läßt sich darum 


— wie alle (auch die nicht-periodischen) Lösungen — in konvergente Potenzreihen nach A ent- 
wickeln. Das Gegenteil kann nur stattfinden, wenn sich keine endlichen Randbedingungen 
(3.5) ergeben, indem die Frequenz (3.10) den Wert Null hat. Für >/* wird also in (3.10) 
©>0 gehen, und dann hört die Konvergenz der Reihe (3.9) auf. 


$ 4. Die erste Näherungsstufe bei fehlender Grenzresonanz. Wenn man die Reihen (3.9), 
(3.10), (3.11) in die Differentialgleichung (3.8) einsetzt, so erhält man in der Grenze M>0 


Das ist die Gleichung des linearen Grenzfalls; ihre allgemeine Lösung ist nach (2.4) 


mit den durch (2.7) bestimmten Grenzfrequenzen @,,, ®g. Da nach (3.9) für M>O die 
Lösung z in z, und z, in die Gestalt (4.2) übergeht, so muß nach (3.5) die durch (4.2) erklärte 
Funktion periodisch in t mit der Frequenz &, sein. &, bedeutet dabei eine bestimmte — an 
sich beliebige, im folgenden jedoch festzuhaltende — der beiden Grenzfrequenzen @g1, @ss- 

Es werde nun vorausgesetzt, daß die andere der beiden Grenzfrequenzen nicht ein un- 
geradzahliges Vielfaches der gewählten Grenzfrequenz &, sei; wenn also w, die eine aus- 
gesuchte Wurzel der Gl. (2.6) ist, so soll keine der anderen Wurzeln die Form (?r +1) o, 
haben: 

+! +’ +a (n=1233,..) . .. (43). 

Um einen kurzen Ausdruck zu haben: es soll keine „Grenzresonanz“, d. h. (ungeradzahlige) 
Resonanz zwischen den Grenzfrequenzen, stattfinden. 

Unter dieser Voraussetzung (4.3) muß, damit (4.2) die Frequenz &, habe, der Amplituden- 
faktor der von dem ausgewählten &, verschiedenen Grenzfrequenz verschwinden; und die 
anderen in (3.5) enthaltenen Forderungen führen dazu, daß in der Grenze M>0 


mut C,=1 


werden muß. Nun ist 2, durch diese Forderung nur bis auf Glieder der Größenordnung A 
bestimmt ''), wie überhaupt alle 2, der Reihe (3.9). Diese Unbestimmtheit läßt sich dazu aus- 
nutzen, den Periodizitäts- und Symmetriebedingungen, die durch (3.5) der ganzen Lösung z 
auferlegt werden, auch durch die einzelnen Koeffizientenfunktionen 2, zu genügen, so daß ihr 
Erfülltsein durch die Gesamtfunktion sich daraus als selbstverständliche Folge ergibt. Es ist 


dann zu setzen 


Dabei ist von der Frequenz » bisher allerdings nur das erste Glied der Reihe (3.10) bekannt. 
Für 2,,2,,... sind also jetzt die Randbedingungen: 
Die Glieder bis zur Größenordnung 4' in (3.8) geben die Gleichung 
+0), +a +4, +(a+c)2, + a’ 


4.6). 
en) =0 (4.6) 

Nun folgt aus (4.4): 

45 
Aus (4.6) wird also 
— (a? + c”) + a? b?) cos w t 
1 

0083 + 3 cos wt — a? b? cos 3. — 3a? b? cos wi) 


11) Handbuch der Physik, Bd. V. Berlin 1927, S. 154. 
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Die zu den Randbedingungen (4.5) gehörige Lösung von (4.7) muß offenbar die Form haben 


(49). 


Führt man (4.8) und die Reihe (3.10) in (4.7) ein, so geben die Glieder mit dem Faktor 4! (die 
mit 4° verschwinden zufolge (2.6)) die beiden Gleichungen 


mit 


403-2(@ + +2 a W)+0:0,=0. 


Der Faktor von C,, in (4.10.1) ist die linke Seite von (2.6), also null; der Faktor von (,, in 
(4.10.2) hingegen verschwindet zufolge der Voraussetzung (4.3) nicht, so daß C,, sich aus 
(4.10.2) bestimmen läßt. C,, folgt dann aus (4.9); die Gl. (4.10.1) liefert dafür den Frequenz- 
koeffizienten 


(a? + 3.08) 160, (4.11), 
wi; —a?b? 


Wenn e?>a?’+b?, was praktisch oft vorkommt (hohe Drehzahl des Kreisels = starke Kopplung), 
so lassen die Gl. (2.7), (4.11), (4.12) sich in erster Näherung ersetzen durch 5 


Präzession Nutation 

b /a\? 

1 1 


Das heißt: Die Präzession ist in der Frequenz, die Nutation in der Schwingungsform ziemlich 
unabhängig von der Schwingungsweite M. Dagegen ändert sich die Form der Präzession in 
dem Sinne, daß an. Stelle der harmonischen Schwingung mehr eine gleichförmige Bewegung 
mit ruckhaftem Richtungswechsel tritt, an Stelle der Sinuslinie eine Zickzackkurve. Die Frequenz, 
besonders die der Nutation, wird — wie zu erwarten war — mit wachsendem Ausschlag kleiner. 


$ 5. Die weiteren Näherungen bei fehlender Grenzresonanz. Die Glieder der Größen- 
ordnung #? führen wegen der in (3.11) auftretenden Produkte 


10 
für 2, zu dem Ansatz 
. 2.2... BU). 
Eine nähere Betrachtung der Reihen (3.7) bzw. (3.11) lelırt, daß es so weitergeht: 2, wird aus 


vier Kosinusgliedern der Frequenzen ®, 3®, 5o, 7 bestehen, in z, tritt außerdem noch die 
Frequenz 9 auf, usf. Man kann daher allgemein ansetzen 


n 
wobei zur Befriedigung von (3.5) noch 
n 
0,,=0 


gelten muß. — Schreibt man die Reihe (3.11) als 
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so enthält S,„(z) immer den Summanden 2, und außerdem nur Produkte der 2, mit m<n; 
und zwar so, daß auch S„(z) als eine Kosinussumme, mit denselben Frequenzen wie in z,, 


darstellbar ist: 


In dieser Formel ist nach dem eben Gesagten der Koeffizient D,„, eine Summe aus dem 
Koeffizienten C,„, der Reihe (5.2) und Produkten der C„. mit m <n; diese weiteren Glieder 
sollen mit E„, bezeichnet werden: 


Die Ausdrücke E,„, sind in der Tafel 1 für n=1,2,3,4 aufgeführt '”): 


Tafel 1. Tafel der E„, in der Gl. (5.7). 


Eu 
1 
Eu = s 
1 
1 
En— 
1 1 
Eun=+ g 
1 


Ct 


RER! 3 5 1 


1 3 1 1 
1 Et 1 1 1 
1 
3, Br 3 1 1 1 1 1 
Eu=— g 53 040 Cw + 737 380 
1 1 1 
8 8 4 (Ca 3 38,0 +6C,+4C%) 1 105 920 
1 1 1 


Die periodischen Lösungen (3.9) lassen sich in jedem ganz im Innern ihres Konvergenz- 
gebietes gelegenen abgeschlossenen Bereich durch ihre Fourierreihen darstellen; diese Fourier- 
reihen erhält man durch formale Umordnung der Doppelreihen (3.9), (5.2): 


‚=® 


12) Ein Weg zur Berechnung der E„» ohne Benutzung der S„(z) ist im $ 11 angegeben. 
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Ebenso kann man S(2) darstellen durch 


n=r 


der Zerlegung (5.7) entspricht hier die Zerlegung 


Werden die Reihen (5.8), (5.10) in (3.5) eingesetzt, so findet man zur Bestimmung des 
Koeffizienten C, von cos (2»+1)ot die Gleichung 


Darin sind zur Abkürzung geschrieben worden 


Alle vier Größen in (5.14) sind analytische Funktionen von 4. Die Reihenentwicklung von 
A, und 42,, die man durch Benutzung von (3.10) erreicht, sei 


. (5.18). 
| 
Für die Koeffizienten A„, in (5.17) bieten sich einfache Ausdrücke: 
. 5.19). 
Av=@r m | 


Gleich einfache Formeln für die Koeffizienten 2,„, sind wegen der in (5.16) stehenden vierten 
Potenz der Reihe (3.10) für & nicht möglich. Doch wenn man setzt 


so gelten für die ersten » nachstehende Darstellungen: 

+1? 

2, = +1 [2 + 

2, =2r+1)[2w, ©, +2 + 

2, = +1)’ [2 ww, + @, + 

2, = +1’ Po, +20, + @}) + 
(120,0, + 

(0,0, + +40} 


Da man die Reihenentwicklung (3.9), (3.10) irgendwo abbrechen wird, so hat man dann 
einen Reihenabschnitt, der eine Näherung für die vollständige Lösung darstellt; und die 
Hinzunahme der Glieder mit der nächsthöheren Potenz von A ist eine neue Stufe auf diesem 
Näherungswege. Auf der n-ten Näherungsstufe handelt es sich um die Bestimmung von @,, 


. 6.21). 
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Cyor Onis +: Cnn- Den Koeffizienten C„, berechnet man, indem man in (5.14) die rechte 
Seite nach A entwickelt und den Faktor von 4” nullsetzt: 


n—v 
(weiter als bis na — » braucht nicht summiert zu werden, da dann die Faktoren E und (€ ver- 
schwinden). In A„, und in 2, treten nur die Frequenzkoeffizienten &, bis &„ auf, in den 
En, nur die C„,„ mit h<m. Für »>0 enthält die Gl. (5.22) mithin außer dem gesuchten 
C,„, nur Größen, die schon bei den voraufgehenden Näherungen ermittelt wurden. Der Faktor 
von ('„, ist Q,,; nach der eingangs gemachten Voraussetzung (4.3) ist aber 


So fließt aus (5.22) die Endformel für die C,, # >0): 


C„. kann dagegen nicht aus der Gl. (5.22) für »=0 entnommen werden, weil sein Faktor 
2,, zufolge (2.6) verschwindet; vielmehr ist hierfür die Gl. (5.3) heranzuziehen: 


n 
Statt dessen liefert die erste Gl. (5.22) den neuen Frequenzkoeffizienten w,„, der dort in 2,, 
in Gestalt des Summanden 2, @,„ 4, vorkommt. Alle anderen Größen sind bereits bekannt; 
setzt man 


so hat man unter Beachtung von E,,=0, C,,=1: 
n—]1 n—1 


Unter der Voraussetzung (4.3) fehlender Grenzresonanz also liefern die Gl. (5.2), (5.24), (5.25), 
(5.27) stufenweise alle Glieder der Reihen (3.9), (3.10), die die gesuchten periodischen Lösungen 
darstellen. 


$ 6. Die periodischen Lösungen bei Grenzresonanz. Nur wenn &, die Grenzfrequenz der 
Präzessionen bedeutet, kann Grenzresonanz eintreten — für die Nutation ist immer die in 
$ 4 und $ 5 angegebene Entwicklung möglich. Es seien jetzt die Bedingungen (4.3) = (5.23) 
nicht erfüllt, sondern es sei 


Nach (5.21) ist dann 


und nach (2.7.1) folgt daraus 


die „kritische“ Kopplung ky hat somit die Größe 

Diese Wurzel ist nur reell, wenn 


Es gibt bei festen Werten von a? und b? eine abzählbare Menge von kritischen '*) Kopplungs- 


13) „Kritisch“ nicht eigentlich in physikalischem, sondern nur in mathematischem Sinne, sofern die Unter- 
scheidung zwischen reinperiodischen und fastperiodischen Lösungen weniger physikalisch als mathematisch ist — s. u.$8. 
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zahlen ky; der Index N, der kleinsten kritischen Kopplung ist durch die Ungleichung (6.3) 
bestimmt. 

Als Beispiel für den Gang der Lösung unter diesen Umständen möge zunächst der ein- 
fachste Fall N=1 behandelt werden. Da Q,,=0, so versagt (5.24) für alle C,„, — aber schon 
die Formel (4,4) für 2, wird ungültig. An ihre Stelle tritt jetzt als die allgemeinste in (4.2) 
enthaltene, den Bedingungen (3.5) genügende Formel 


mit der Bedingung 


wobei wie oben im $ 4 nicht nur z, sondern auch 2, die Periodizitäts- und Symmetrie- 
forderungen (3.5) erfüllt. Das Verhältnis von C,, und C,, bleibt aber auf dieser nullten 
Näherungsstufe noch unbestimmt. 


Auf der ersten Stufe erhält man für S,(z) jetzt 


(6:6) 


mit folgenden Werten der E,,: 


1 


1 


2, wird nun auch die Frequenzen », 30, 5o, enthalten: 
2, . . (68). 


Allgemein ist der Unterschied gegen den Fall ohne Grenzresonanz der, daß die Funktionen z, 
mehr Teilschwingungen enthalten. Die Reihen (5.8), (5.10) bleiben bestehen; in (5.2) ist formal 
bis = zu summieren (in Wirklichkeit treten bei jedem » nur endlichviele Glieder auf, 
deren Anzahl jedoch auch noch von N abhängt). In den Gl. (5.9), (5.13) ist von n=0 an zu 
summieren (wobei in Wirklichkeit i. A. viele Anfangsglieder ausfallen). Die Gl. (5.14) bleibt 
in Kraft, doch in der daraus gewonnenen (5.22) muß m jetzt von O0 bis n laufen. Mit der- 
selben Änderung bleibt für vz&0, N auch die Formel (5.24) in Geltung; an der Gl. (5.27) 
ändert sich nichts. In (5.3) ist wie in (5.2) bis v„=»x zu summieren. 

C,, und C,, sind jetzt bestimmbar. Die Gl. (5.22) lautet für »=0,1 (d.h. die beiden 
Werte, wo (5.24) versagt): 


In beiden Gleichungen tritt außer C,,, C,, — zwischen denen die Gl. (6.5) besteht — noch 
der Frequenzkoeffizient &, auf, und zwar nach (5.21) in (6.9.1) mit dem Faktor 2, 4,, in 
(6.9.2) mit dem Faktor 180, 4A,,. Aus der Voraussetzung (6.1) mit dem daraus abgeleiteten 
Werte von ®; folgt nun aber für N=1, daß 4, ie und allgemein 


Wenn man also (6.9.1) durch C,, kürzt — auch E,, enthält nach (6.7) diesen Faktor'*) — und 
nach Multiplikation mit 9 C,, zu (6.9.2) addiert, so hebt sich die unbekannte Frequenz w, 
heraus; übrig bleibt eine Gleichung 3. Grades für C,, und C,,, aus der sich diese zwei 


14) Bei beliebigem N ebenfalls: denn Coo=0 (und auf den weiteren Stufen Cyo=0) ist immer eine mögliche 
Lösung von (6.9), und zwar ist es die hier nicht interessierende mit > wg. 
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Koeffizienten unter Beachtung von (6.5) berechnen lassen. Da die Gleichung vom 3. Grade 
ist, so hat sie stets eine reelle Lösung. Durch (6.7) sind jetzt die E,, bekannt und damit 
dank den Formeln (5.24), (5.27) auch C,,, C,,, C,, und o,. 

Unbekannt bleiben auf dieser Näherungsstufe noch C,, und C,,; pur ihre Summe ist 
durch (5.3) schon sicher. Diese beiden Koeffizienten werden nun in gleicher Weise dazu ge- 
braucht, um auf der zweiten Näherungsstufe die beiden Gl. (5.22) für n=2, v=0,1 
zu befriedigen, in denen wiederum C,, und C,, nur mit verschwindenden Koeffizienten auf- 
treten. &, kommt in diesen beiden Gleichungen nur linear vor, läßt sich also wieder eliminieren. 
Es bleibt eine Gleichung 5. Grades mit mindestens einer reellen Lösung. Sind C,, und (C,, 
gefunden, so folgen C,,, ... C,, auf dem gewöhnlichen Wege; usf. 

Entsprechend ist zu verfahren, wenn N >1 ist. Für z, nimmt man den Ausdruck 
N+1)wt mit vorläufig unbestimmten Koeffizienten C,,, C,y, deren 
Summe 1 sein muß. z, wird ein Kosinuspolynom aller in S, (2) — 2, auftretenden Frequenzen. 
Auf der n-ten Stufe kommen die Koeffizienten C„, und C„y nur mit den verschwindenden 
Faktoren (,, und @,, vor; dafür lassen sich aus diesen beiden Gleichungen, nachdem der 
in ihnen linear auftretende neue Frequenzkoeffizient &,„ eliminiert worden ist, die beiden auf 
der vorigen Stufe unbestimmt gebliebenen Koeffizienten C,„_,,, und C„_,,y durch eine Gleichung 
ungeraden Grades berechnen, da ihre Summe schon bei der (n — 1)-ten Näherung ermittelt 
wurde. Sind die Werte dieser beiden Koeffizienten festgestellt, so ergeben sich auch &„ und 
die C„, mit Ausnahme der auf der nächsten Stufe erst zu bestimmenden C,„, und C,„x- 

Damit ist der Verlauf der Berechnung der periodischen Lösungen im Falle der Grenz- 
resonanz umrissen. 


$ 7. Die periodischen Lösungen des Systems (2.2). Die Gl. (3.8) für z ist durch das Ver- 
fahren der $$ 3bis6 gelöst. z hat die Gestalt einer Potenzreihe (3.9), deren Koeffizienten z, 
Kosinuspolynome (5.2) sind, oder umgeordnet die Gestalt einer Kosinusreihe (5.8), deren 
Koeffizienten Potenzreihen (5.9) nach dem Parameter 4 sind. Die Grundfrequenz & ist gleich- 
falls durch eine Potenzreihe (3.10) gegeben. Dieses 2 löst zugleich das System (3.4) bei den 
verlangten Randbedingungen (3.5). Um die andere in (3.4) eingehende Funktion « (f) zu be- 
stimmen, eliminiere man aus (3.4) die Funktion S (2); es bleibt 


Man sieht durch Einsetzen der Reihe (5.8) für 2, daß « die Form einer Sinusreihe haben muß: 


ihre Koeffizienten ergeben sich unter Benutzung der Bezeichnung (5.15) als'*) 


Durch Multiplikation mit der Konstanten M (3.2) findet man dann die periodischen Lösungen 
x(t), y(t) des Systems (2.2) zu den Randbedingungen (3.1). Macht man die Substitution (2.1) 
rückgängig, d. h. multipliziert man & (2) noch mit Y, so ist auch das ursprüngliche System (1.3) 
gelöst; der Rollwinkel y des Schiffes hat also die Sinuskoeffizienten M YB,, der Pendel- 
winkel # des Kreisels hat die Kosinuskoeffizienten M-C,. Doch soll den weiteren Betrach- 
tungen das System (2.2) zugrundeliegen. 

Es sind jetzt zwei Scharen periodischer Lösungen von (2.2) bekannt, die zu den Grenz- 
frequenzen und gehören: 


Diese Lösungen und ihre Frequenzen hängen ab von der Schwingungsweite M der Funk- 
tion y(t). Da in (2.2) die Zeit nicht explizit auftritt, so können die Lösungssysteme von (2.2) 
beliebig in der Zeit verschoben werden. 

Ebenso wie beim linearen System (2.3) stellen sich auch hier die periodischen Lösungen 
nur bei passenden Anfangsbedingungen ein, sonst ist die Lösung i. A. nicht periodisch '®). 


15) Aus (7.3) ließe sich leicht die Potenzreihe By = Ni" B„» ableiten, falls dies gewürscht wird. 


16) Poincare hat a. a. 0.9) bewiesen, daß es auch reinperiodische Lösungen gibt, die wohl der Form nach 
für > in die periodische Lösung des linearen Problems übergehen, aber nicht in der Frequenz; statt dessen gilt 
o> ” (m=2,3,4,..). In nullter Näherung tritt in diesen Lösungen also nur die m-te Harmonische der Grundsehwingung 
auf, während die Amplitude der Grundschwingung selbst den Faktor X enthält. Darum haben diese Lösungen weniger 
Interesse, so daß ihre Ausschaltung durch die Forderung (3.1) >, berechtigt war. 
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$ 8. Die nichtperiodischen Lösungen. Das Gleichungssystem (2.2) in der Form (3.6) 
gehört'”) zu den Lieschen Systemen mit Fundamentallösungen **). Als Fundamentallösungen 
kann man z. B. die oben ermittelten reinperiodischen Lösungen wählen. Dann erkennt man, 
daß die allgemeine Lösung zu den „im weiteren Sinne periodischen“ Bohlschen Funktionen 
gehört'*), d. h. zu den fastperiodischen Funktionen, deren Modul eine Basis aus endlichvielen 
Frequenzen hat. Diese Basis muß mindestens zwei Frequenzen enthalten, da sie stetig in die 
der Lösung (2.4) des linearen Systems übergeht. Wendet man die Methode der schrittweisen 
Näherungen auf (3.6) an, indem man von der linearen Lösung 


y=A,coso, t+B, sino, t+ A,cos B,sin m, 


ausgeht, so entsteht für sin y eine Fourierreihe der Form 


sin y= (Acoswt+ Bsin wi), 


wobei alle Werte 


durchläuft. Dieselbe Frequenzenmenge tritt darum in der zweiten Näherung für y auf. In den 
folgenden Näherungen treten keine neuen Frequenzen mehr auf, da die Fourierzerlegung des 
Sinus einer trigonometrischen Reihe mit den Frequenzen (8.1) nur wieder diese Frequenzen 
enthält?°). Das System (2.2) läßt sich also formal durch trigonometrische Reihen dieser Art 
lösen; und da die allgemeine Lösung zu den Bohlschen Funktionen gehört, so muß sie auch 
wirklich diese Form haben: der Modul der Lösung hat eine Basis von zwei Frequenzen, und 
es tritt nur die Hälfte aller Frequenzen des Moduls auf. 

Von hier aus läßt sich auch der Unterschied zwischen den periodischen Lösungen mit 
und ohne Grenzresonanz überblicken. Bei Grenzresonanz gibt es keine reinperiodische Lösung 
von (2.2), die für M>0 in die reine Präzessionsschwingung überginge; sondern eine solche, 
welche in die mit einer Nutationskomponente überlagerte Präzessionsschwingung übergeht — 
diese Grenzform ist dabei freilich auch eine reinperiodische Funktion. Andrerseits gibt es 
bei Grenzresonanz sicher fastperiodische Lösungen, die in die reine Präzessionsschwingung 
übergehen. Denn die Funktion, nach der sich die Fundamentallösungen zur allgemeinen 
Lösung zusammensetzen, muß aus Stetigkeitsgründen bei kleinen Ausschlägen (der Einzel- 
schwingung wie der zusammengesetzten Schwingung) in die lineare Zusammensetzung übergehen. 


$ 9. Ein Zahlenbeispiel. Zur Veranschaulichung der Konvergenz des oben entwickelten 
Verfahrens zur Berechnung der periodischen Lösungen von (1.3) bzw. (2.2) wurde ein von 
Klein und Sommerfeld°') angegebenes Beispiel durchgerechnet. Danach waren 


H=0.4m, J — 550000 mkgsek?, 
N=f-T, T=175mkgsek’, f=189sek". 


Daraus folgt für die in (2.1) eingeführten Konstanten: 
Y—0.016514, a—=0.6182, k’=13.2601. 


Da h nicht bekannt ist, ließ sich ein Wert für die ungekoppelte Schwingungsfrequenz b des 
Kreisels nicht angeben. Doch wird a. a. O. S. 841 die Größe v’—b’/a? als wahrscheinlich 
etwas kleiner als 1 bezeichnet. Der Rechnung hier seien darum folgende abgerundete Werte 
zugrunde gelegt: 


V=00165, «=062, k’= 13.26. 
Dieses k? ist nicht kritisch; die nächstgelegenen kritischen Werte der Kopplung sind nach (6.2): 
k,=123%8 und 13.574. 


17) Dort ist dann noch M=4=1 zu setzen: so wird S(%,)=sin y, =X. 

18) Enzyklopädie der math. Wiss. Bd. IIla, Leipzig 1899 bis 1916, S. 275 bis 276; S. Lie: Über Differential- 
gleichungen, die Fundamentalintegrale besitzen. Ber. üb. d. Verh. d. Kgl. Sächs. Ges. d. Wiss. zu Leipzig, Math. 
Phys. Kl., Bd. 45 (1893), S. 341 bis 348. 

19) Piers Bohl: Über die Darstellung von Funetionen einer Variabeln durch trigonometrische Reihen mit 
mehreren einer Variabeln proportionalen Argumenten. Dorpat 1893. — Harald Bohr: Fastperiodische Funktionen 
Ergebnisse der Math. Bd. 1 (1932), Heft 5, S. 3. 

20), Der Beweis folgt aus den unten mitgeteilten Formeln (11.2). 

21) A. a. 0.1), S. 532. 
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Die Ordnung (11 bzw. 12) dieser kritischen Zahlen ist eroß genug, um eine gute Darstellung 
der Lösung durch die ersten Glieder erwarten zu lassen, obwohl ein monotones Abnehmen 
der Reihenglieder (dem Betrage nach) vielleicht erst jenseits von n=12 eintreten wird. 
Auf 6 Stellen abgerundet sind die Werte der C,„, für n=0,1,2,3 in der Tafel 2 an- 
gegeben, dazu die ®„. Die eingeklammerten Angaben der ersten beiden Stufen sind die mit 
den Näherungsformeln (4.13) erhaltenen Werte. Die Näherungsformeln sind bei dem vor- 
liegenden Verhältnis e?:(a?+ 5°) — 11,7 schon gut verwendbar. Da bei der Nutation die 
Frequenzreihe ziemlich langsam konvergierte, ist dort auch noch ®, aufgeführt worden. — 


Tafel 3. Frequenz und Einzelamplituden der periodischen 
Lösungen des Systems (1.3), (2.2) für die Ausschlaggröße 
M = 60°: 


Tafel 2. Tafel der Entwicklungskoeffi- 
zienten für die Frequenz und die Teil- 
schwingungsweiten 


der Präzession: der Nutation: Nutation 


Präzession 


T 
| ?—=354, sec! 


% = 0.155192 3.796 830 0.15477 see”! 
(0.158 502) (3.717 526) T® = 40.597 sec 
C,=0946 +56.8° C,= + 1.0052 [+ 60.31° 
©, = —0.0W402 — — 0.227 476 C,=0.047 28 C, = — 0.0053 / — 0.32 
(— 0.000 444) (— 0.232 345) C,=0.007 I+ 0.4 C, + 0.0000 / 0.00° 
00 — 0.0000 / — 0.00% 
— 0.041 667) .005 2 _— — 
Cu= — — 0.01447 | — 0.87° Bo + 0.01550 / + 0.93° 
PB,=+ 0.00011 /+0.01° PB, = — 0.00024 / — 0.01° 
— + 0.00002 / +0.00° B, — 0.00000 / + 0.000 
— + 0.00000 / + 0.00° — 0.00000 / — 0.000 
+ 0.004 515 + 0.000 046 Frequenz und Amplituden der periodischen Lösungen des 
Systems (1.2), (2.3) für dieselbe Ausschlaggröße M = 60°: 
= +0.0004 — — 0.001 7236 
Co= — 0.001307 + 0.000 022 0.15519 sec”! 3.179683 sec”! 
— 0.0033 — 0.000 018 = 40.487 sec Ta = 1.655 sec 
+000095 — 0.000 004 
Ca = + 0.000 695 Oz; = — 0.000 000 C,= 1.000 / + 60.0° C,= 1.0000 / + 60.00° 
— 0.000 402 PB, — — 0.01565/— 0.9 | + 0.01654 — 0.99° 


Aus diesen Koeffizienten lassen sich Frequenz » und Teilamplituden C, für verschiedene 
Ausschläge M berechnen. In Tafel 3 sind die zu M=60° = 1.047 gehörigen Werte zusammen- 
gestellt. Die Gradangaben sind die entsprechenden Winkelamplituden der in den Funktionen 
d(t) und y(t) des Ausgangssystems (1.3) enthaltenen Teilschwingungen. T=27[/o ist die 
Schwingungsdauer in Sekunden. Zum Vergleich mit (1.2) sind unten noch die zum gleichen 
M gehörigen Zahlenwerte aufgeführt, die sich aus dem linearen Ansatz (2.3) ergaben. 

Führt man diesen Vergleich durch, so zeigt sich: Die Frequenz der Präzession hat um 
0.3 vH. abgenommen, die der Nutation um 6.6 vH. Diese Frequenzen würde das lineare 
System (2.3) liefern, wenn dort bei gleichbleibendem a die beiden Konstanten b und k um je 
6.9 vH. vermindert würden. Für (1.3) heißt das, daß die „effektive“ Umdrehungszahl f des 
Kreisels um 6.9 vH. geringer ist als die tatsächliche Umdrehungszahl und daß seine Schwingungs- 
höhe h effektiv um 13.3 vH. kleiner ist als bei sehr kleinen Ausschlägen. 

Als Maß der Formänderung könnte man die Änderung des Koeffizienten C, der Grund- 
schwingung ansprechen. Diese Formänderung beträgt bei der Präzession — 5.4 vH., bei der 
Nutation +0.5 vH. Diese Formänderung bedeutet, daß die Beschleunigung der Kreisel- 
schwingung im Umkehrpunkt (£=0), für welche sich aus (5.8) mit (5.9) und (5.3) die Formel 


+90, +350,+49 0, 


ergibt, gegenüber einer reinen Sinusschwingung von gleicher Frequenz und gleichem Ausschlag 
bei der Präzession auf das 1.64-fache zugenommen hat; dem steht bei der Nutation eine Ab- 
nahme um etwa 4 vH. gegenüber. Bei Berücksichtigung der Frequenzänderung kommt zu 
diesen Angaben noch der Faktor ®’/w}; demnach ist die Beschleunigung im Umkehrpunkt 
bei der Präzession um 62 vH. größer, bei der Nutation um 16.4 vH. geringer, als sie es naclı 
dem linearen Gesetz (1.2) bei gleicher Schwingungsweite M=60° wäre. 


(9.1) 
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Die zu diesem Kreiselausschlag von 60° gehörigen Schlingerwinkel des Schiffes, die man 

aus (2.1), (3.2) und (7.2) für » I=4 2 zu 

0.92) 
findet, sind bei der Präzession um 7.0 vH., bei der Nutation um 4.8 vH. kleiner als nach der 
linearen Theorie. Zu gleichen Schiffsausschlägen gehören deshalb größere Kreiselausschläge, 
als es linear der Fall wäre. 

Insgesamt ist demnach festzustellen: Soweit es auf die Eigenschwingungen des Systems 
Kreisel-Schiff unter Außerachtlassung der Dämpfung ankommt, ist bei Kreiselausschlägen bis 
zu 60° (die etwa die Grenze der wirklich vorkommenden angeben) eine größere Abweichung 
der linearen Gleichungen (1.2) von den nichtlinearen (1.3) nicht festzustellen — die Diffe- 
renzen liegen innerhalb von 10 vH. Lediglich die Beschleunigung der Präzession in den 
Umkehrpunkten wird nahezu um die Hälfte zu klein geliefert (die der Nutation um ein Fünftel 
zu groß). Die Beanspruchung der Schiffsverbände durch den Umkehr-Ruck ist also doppelt 
so groß wie nach der linearen Theorie. 

Die oben berechneten Näherungsstufen reichen noch nicht aus, um die Verhältnisse bei 
einem größeren Ausschlagwinkel, etwa M=%", mit derselben Genauigkeit darzulegen. Dazu 
wären bei den Amplituden noch etwa zwei Stufen notwendig. Die Frequenzen nehmen die 
Werte &" = 0.1542, &® —3.2%0 an, werden also um 0.6 vH. bzw. um 15.8 vH. kleiner als die 
Grenzfrequenzen. Das entspräche bei Weiterverwendung des linearen Ansatzes (1.2) einer 
Verringerung der wirklichen Umdrehungszahl f des Kreisels um 16.5 vH. und seiner Pendel- 
länge h um 29.9 vH. — Über die Vergrößerung der Umkehrbeschleunigung bei der Präzession 
sind an Hand der bisher ausgerechneten Koeffizienten auch rohe Angaben noch nicht möglich. 


$ 10. Das ursprüngliche System (1.1). Die periodischen Lösungen des Systems (1.1) 
können grundsätzlich in gleicher Weise gewonnen werden. Nach Division mit J, j bekommt 
man an Stelle von (2.2) die Gleichungen 

. (10.1). 
j+b’siny+ksinz—0 


Sie unterscheiden sich von den Gl. (2.2) in folgendem: Es besteht nicht die Möglichkeit, 
durch eine lineare Transformation wie in (2.1) die beiden Kopplungskoeffizienten gleichzu- 
machen. Außerdem läßt sich keine der beiden Veränderlichen «, y eliminieren, ohne daß 
eine unbequeme transzendente Funktion aufträte. Deshalb ist gemeinsame Weiterbehandlung 
beider Gl. (10.1) vorzuziehen; das führt dazu, daß, wo früher eine lineare Gleichung auf- 
zulösen war, jetzt immer ein System von zwei linearen Gleichungen zu lösen ist. — Weiter 
ist es wohl unvermeidlich, daß die Gleichberechtigung von x, y in (10.1) durch die Wahl 
des die Schwingungsweite kennzeichnenden Parameters M gestört wird. Man kann z. B. 
wieder den maximalen?) Ausschlag der Veränderlichen y mit M bezeichnen. Die nach dem 
Stetigkeitssatz von Poincare für hinreichend kleine M bestimmt existierenden Lösungen 
sind wieder in absolut und gleichmäßig konvergente Fourierreihen entwickelbar: 


x—=MB,,sinot+ MB, ,sin3ot+ MB,,sin5ot-+.. 
und ebenso die Sinus dieser Funktionen: 
sinz= 
. (10.3). 


siny= MC}, cosot+ MC],„cos3wt-+... 
Der Zeitnullpunkt wurde wieder so gelegt, daß für t=0 y gerade seinen positiven Extrem- 
wert und x den Wert O hat. Die Koeffizienten B,, C,, B;, C, (die letzten sind in (5.10) mit 
D, bezeichnet worden) sind ihrerseits durch Potenzreihen nach A= M? darzustellen, ent- 
sprechend den Gl. (5.9), (5.11). Wenn keine Grenzresonanz stattfindet (4.3), werden sie wieder 
von der Größenordnung 4” sein, so daß 


nv 


. (10.4). 
C,= Ni" G= 


n=v 


22) Besser wäre der Ausdruck „Mittelausschlag“, denn es ist denkbar, daß bei größerem M die Oberschwingungen 
so stark werden, daß M, obwohl Symmetriepunkt der Kurve, doch relativ zur Umgebung ein Minimum darstellt. 
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Bedeutet M den Extremausschlag von ,, dann wird wieder 


sein. Die C,/, unterscheiden sich von den C„, um die oben angegebenen Ausdrücke E,„,»; 
die B,, von den B,„, um ähnlich gebaute Ausdrücke, die von den E„, außer in einigen Vor- 
zeichen (weil es sich um eine Sinus- statt einer Kosinusreihe handelt°®)) noch darin ab- 
weichen, daß die durch (10.5) ermöglichten Vereinfachungen fortfallen. Zur Veranschaulichung 
seien die ersten dieser Beziehungen einander gegenübergestellt: 


Bi = Bir 
C„=C, 
Einsetzen der Reihen (10.2), (10.3) in (10.1) gibt 


Hier sind wieder die Entwicklungen (10.4) einzusetzen. Auf der nullten Stufe hat man dann 
das lineare System 


Bo —hw, — =0 
Co — kw, Bir — =0 
aus dem sich bei Berücksichtigung von (10.5), (10.6) &, und B,, ermitteln lassen. «, hat die 
auch aus (2.6) bekannten Werte?*), B,, ist im wesentlichen das r, aus (2.5). — Auf der 
ersten Stufe entstehen vier Gleichungen: 

+2%,®, Bo ho, Ch ho, B,=0 

+20,0, 0, ko, Bo Bu —bC,=0 

B,=0 

C,=0 


. (10.9). 


Mit Benutzung von (10.6) lassen sich aus den letzten beiden Gl. (10.9) B,, und C,, berechnen. 
Dann dient (10.5) dazu, C,, zu finden, und die ersten beiden Gleichungen in (10.9) liefern die 
beiden noch fehlenden Unbekannten, nämlich B,, und o,. — Die Auflösung der Bestimmungs- 
gleichungen der zweiten und der folgenden Stufen geschieht ebenso: aus je zwei Gleichungen 
können B,„, und C,, festgestellt werden, solange »>0; dann hat man C,, aus (10.5) zu ent- 
nehmen, und die beiden ersten Gleichungen ergeben darauf B,„, und &„. Allgemeine Formeln 
dafür, die den oben abgeleiteten Formeln (5.24), (5.27) entsprächen, ließen sich ohne Schwierig- 
keit aufstellen. Auch das Lösungsverfahren des $ 6 für den Fall der Grenzresonanz ist ohne 
weiteres übertragbar; und die Ausführungen des $ 8 können so gut wie wörtlich über- 
nommen werden. 


$ 11. Anhang: Die Fourierentwicklung des Sinus einer Fourierreihe. Die in Tafel 1 auf- 
geführten Koeffizienten E„, der 4-Entwicklung stellen diejenigen Größen der Formeln (5.24) 
und (5.27) vor, die die unübersichtlichste Zusammensetzung haben. Nach ihrer Definition 
(5.7) sind sie durch die D„, erklärt; und diese wurden oben (5.5) durch Kosinuszerlegung 
von S„(z) gewonnen, wo S„(z) der Koeffizient von 4" in der Reihendarstellung von S (2) 
war, die sich aus (3.7.1) durch Einsetzen von (3.9) und Umordnung nach den Potenzen von 
i ergab, nämlich in (3.11.1). — Eine andere Möglichkeit zu ihrer Berechnung besteht nach 
(5.11) in der Potenzentwicklung von D,, wodurch die Benutzung der S,„(z) vermieden wird. 
Nach (5.10), (5.8), (3.3) sind die D, Funktionen der C,, und zwar sind M D, die Kosinus- 
koeffizienten von sin y, während MC, die von y sind: 


23) Infolgedessen ist hier (11.9) ohne die Vorzeiehentransformation (11.4) brauchbar. 
24) Dabei ist darauf zu achten, daß nieht das k in (10.1) mit dem in (2.2) identisch ist, sondern daß das Produkt hk 
aus 10.1) gleich dem k? aus (2.2) ist. 
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y=Mz= MC,cos @v+1)wt 


MC,= Mert:c,, 


siny=MS(e2)= (2v+1)wt 


MD= 


Die Bestimmung der E„,, und das heißt zunächst: der D,, läuft also hinaus auf die Fourier- 
entwicklung von sin y, wo y selbst durch eine Fourierreihe gegeben ist (11.1), wenn auch 
eine spezieller Art. Nun gelten bekanntlich?®) die Entwicklungen 


sin (asin)—=2 N (a) sin(2n +1)t 


. (11.2); 
cos (asin )—J,(a) +2 (a) cos t 


darin sind die J„ die Besselschen Zylinderfunktionen der n-ten Ordnung: 


a n+em 


Die beiden Kosinusreihen in (11.1) werden durch 


in Sinusreihen umgeformt: 


Der Koeffizient d,, der Reihe (11.6) ist durch Multiplikation mit sin (2m + 1) und Integration 
nach r von O bis 27 nebst folgender Division mit x zu ermitteln. Nun läßt sich der Faktor 
sin(2m-+1)r aus (11.5) in der Form darstellen 


sin 
benutzt man dies und vertauscht noch die Reihenfolge?®) von Differentiation und Integration, 


so folgt 


in 


0 


Für diese Größe Y läßt sich aus (11.2) bei Beachtung der Additionstheoreme von sin und cos 
die Produktsumme ableiten 


(n) 


+l)n,=0). 


25) N. Nielsen: Handbuch der Theorie der Zylinderfunktionen. Leipzig 1904, S. 65. — Dagegen konnte Verf. 
die unten daraus abgeleitete Formel (11.9) sonst noch nicht finden. Vgl. dazu auch G. Hamel: Über erzwungene 
Schwingungen bei endlichen Amplituden, Math. Ann. Bd. 86 (1922), Seite 1 bis 13, bs, S. 5. 

26) Das ist erlaubt, da die Reihe (11.1) für y als absolut und gleichmäßig konvergent vorausgesetzt werden kann 
und daraus dasselbe für alle anderen Reihen und Produkte dieses Paragraphen folgt. 


16 


AUM 
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Daraus folgt dann für d,„ aus (11.7) mit der Beziehung 
(2) = — In+ 1 
der endgültige Ausdruck: 


(n) 


((Sv<o; 


Das ist zusammen mit (11.4) der gesuchte Zusammenhang zwischen den D, und den (,. 
Sollen nun etwa E,, E,,, E,, angegeben werden, so braucht man dafür aus d,, d,, d, 
die Glieder mit M°, wie sich aus (11.1) und (11.4) ergibt. Da ce, von der Ordnung M?”*! 
ist, so ist J„(c,) von der Ordnung M'*'-@”*»; für den geforderten Zweck brauchen aus (11.9) 
für die Werte m =0, 1, 2 nur die Glieder mit 


berücksichtigt zu werden. Es kommen daher nur folgende Kombinationen in Betracht: 


| | Daraus sind dann die Vorzeichenverteilungen festzu- 
m m m stellen, die die Summe (2 m +1) ergeben, wie in (11.9) 
1 gefordert: 
2 1 0 m | N, n, N, 
1 0 0 
3 0 0 | +3 0 1) 
2 0 ) 2 0 N) +1 
1 0 0 +2 +1 0 
0 0 +5 0 


Demnach ist für die Bestimmung von E,, der vollständige Ausdruck (11.9) für d, gleichwertig 
(Zeichen: =) der Summe 


2 (So du, ebenso 


Dabei wurde J„,, für J„(c,) geschrieben und ferner benutzt, daß sich J,, von I nur um einen 
Ausdruck der Ordnung M*”** unterscheidet. — In diese kurzen Ausdrücke (11.10) setze man 
jetzt die Anfangsglieder der Potenzreihen (11.3) ein, soweit sie die Größenordnung M° nicht 
übersteigen: 


1 1 1 


Nun sind rechts die c, durch die Potenzreihen (11.1) für die C, zu ersetzen, wobei auf (10.5) 
zu achten ist?”). Greift man aus der so entstehenden Potenzreihenentwicklung die Glieder 
mit der Potenz M° heraus, so folgt sogleich 


1 1 


Gleichzeitig liefern die Summen (11.10) auch alle E„, mit n<2, sie könnten für diesen 
Zweck indessen auf noch einfachere Ausdrücke reduziert werden. 
Bei kleinen » ist die erste Berechnungsweise der E„,, die aus (5.5), sicher einfacher; 
dieser zweite Weg ist aber mindestens zur Kontrolle bei größeren n wegen seiner gänzlichen 
Verschiedenheit sehr geeignet. Zudem läßt sich hier ziemlich leicht einsehen, warum das 
konstante Glied in E„, bis auf den Faktor 2: (— 1)” gleich dem Koeffizienten von «?”*! in 
der Reihe (11.3) für J,,;,(x) ist. 147 


2’) Falls diese Beziehungen gelten. Die angegebene Methode ist unabhängig davon und kann z. B. auch zur 
Berechnung der B},, in (10.4) dienen, s. Anm.?3). 
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Das Hornersche Schema bei komplexen Die Rechnung erfolgt nach dem Schema 
Wurzein algebraischer Gleichungen. 
Bei Stabilitätsuntersuchungen für Systeme mit An-2a a Ad 
mehreren Freiheitsgraden wird man oft auf al- a 
mit komplexen Wurzeln geführt; deren Realteil PAn PAn-ı pa, 
entscheidet nicht nur, ob das System stabil oder 


instabil ist, sondern gibt zugleich ein Maß für 
die Stärke der Stabilität. 


Bei der Wurzelberechnung nach dem Newton- 


schen Verfahren (Berechnung einer besseren 


braucht man an einer gegebenen Stelle x, die 
Werte von g und g’, die gewöhnlich nach dem 
Hornerschen Schema berechnet werden. 


Im folgenden soll eine naheliegende Abart des 
Hornerschen Schemas für komplexe Wurzeln 
beschrieben werden, die vielleicht nicht neu, aber 
zumindest nicht genügend bekannt ist und die 
gegenüber dem ursprünglichen Hornerschen 
Schema!) mit fast der halben Rechenarbeit aus- 
kommt, falls die Koeffizienten der gegebenen 
Gleichung reell sind. 


Vorgelegt sei die Gleichung 


Berechnet werden sollen g (z,) und g’(z,) an der 
komplexen Stelle —=u-iv, wobei und reell 
sind. Dann genügt z, der quadratischen Gleichung 
mit reellen Koeffizienten 


wobei p=2u, g=— (wW + v?) gesetzt ist. Nach 
(2) kann man nacheinander die Potenzen zu", 
reduzieren: 


+4,24 


wobei zur Abkürzung gesetzt ist 


= Ay 
tmP 
a.’ +4,41 944,429 (3). 
+aygq 
Setzt man 


so gilt 

Es ist dann 


1) Fr. A.Willers: Methoden der praktischen Analysis, 
Berlin und Leipzig 1928, S. 205. Dort sind S. 194 bis 212 
auch zahlreiche andere numerische und graphische 
Näherungsmethoden für komplexe Wurzeln beschrieben; 
ferner H. Scheffler, Die Auflösung der algebraischen 
und transzendenten Gleichungen, Braunschweig 1859, S. 99. 


Bei dem ursprünglichen Hornerschen Schema 
für —=u-+iv dagegen 


Ay Au-ı 
A„u +ia,v 
a, 


utia,_,®v utia,'v 


a, + id, 
a; un 
a ayu-ia,v. 
— Br vtißd,u — 


hat man an jeder Stelle die doppelte Anzahl 
von Produkten auszurechnen und außerdem mehr 
Fehlermöglichkeiten durch das Rechnen mit der 
Zahl i. Im neuen Schema rechnet man ganz im 
Reellen und erst bei Anwendung der Gl. (5) geht 
die Zahl i ein. 

Bildet man aus g,(z) ein Polynom 9(z) in 
gleicher Weise, wie man aus g(r) das Polynom 
9, (2) bildete, so ergibt die Division 


— pr—gq +a”c+a,” . (6), 
wobei die a,” die Koeffizienten von 9, (x) sind: 

92 (2) = a,” + a, 
Einsetzen von (6) in (4) und Differentiation liefert 


mit (T). 
b,=2a,"+a,"pund b),=a, +2a,”g—a,”p 
Das Rechenschema, das sich unmittelbar an das 
obige Schema für die a;’ anschließt, lautet jetzt 


q Ay q u” q a ” 
pa 
p 
p 2 q 
b, bu 


Dabei geben die beiden geschweiften Klammern 
an, welche Ausdrücke zur Bildung von b, und db, 
zu summieren sind. 

Bei Polynomen nicht zu hohen Grades wird man 
g’(z) durch gliedweises Differenzieren aus g (x) 
gewinnen und g’(z,) nach (3) (4) berechnen, und 
nur bei Polynomen hohen Grades wird man das 
eben angeführte erweiterte Schema verwenden. 

16* 
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Zahlenbeispiel. 
= + 0,643 2° — 2,619 2 +3 2? — 0,063 2? — 1.101 2 1,012 =0. 
Als Näherungswert wird 
x, — 0,429 -+ 0,695 i 


verwendet. Dann ist 


0,858 2, — 0,66707. 
Nach dem bisher üblichen Hornerschen Schema würde man rechnen 
2, = 0,429 + 0.695 i 


— 0,063 — 1,101 1,012 
0,4897 + 0,7934 5 0,5971 + 0,9674 — 0,2751 — 0,4456 i 
0,9652 — 0,5958 i — 0,1373 -1- 0,0847 i — 0,7312 4 0,4513 i 
1,3919 + 0,1976 — 0,6412 + 1,0521i 0,0057 + 0,0057 


1 0,643 
0,429 + 0,695 i 


— 2,619 3 

0,4599 + 0,7450 i — 1,1335 — 1,8368 i 
— 0,4830 + 0,2982 i — 0.7250 + 0,4475 i 
1.072 1.0.6953 — 2,6421 1,0432i 1,1415 — 1,3888 


0,429 + 0.695 i 0,6439 + 1,0432 — 1.2717 — 2,0601 i — 0,8557 — 1,3863 5 0,9152 + 1,4827 i 
— 0,9661 + 0,5963 — 1,8645 + 1.1509 1,5971 — 0,9859; 1,5114 — 0,9389 i 
1 1,501 + 1,3903 — 2,9643 + 2,6827 — 1,9947 — 2.298051 2,1333 — 217465 1,7854 + 1,6019; 


9 (2) = 0,0057 + 0,0057 i 
g' (&,) = 1,7854 + 1,6019 i. 


Nach dem hier beschriebenen Schema würde man rechnen: 
0,858 2, — 0,66707 


1 0,643 — 2,619 3 — 0,063  — 1,101 1,012 
— 0,6671 — 1,0013 1.3329 — 0,1896 — 1,0098 
0,858 1,2879 — 1,7144 0,2439 1,2989 
1 1,501 —19082 "0.2843 15138 0,083 0,0022 = 0,0083 x, + 0,0022 
— 0,6671 — 1,5736 0,4278 0,0057 + 0,0057 i 
0,858 2,0240 — 0,5502 
ı 2359 — 0,6413 1,8395 1,9416 
1,9416 — 1,6659 
— 1,5783 2,4542 = 2,3049 x, + 0,7966 
2.3049 0,7966 — 1,7854 -+ 1,60193. 


Die Berechnung von g’(x,) hätte auch durch Bildung der Ableitung erfolgen können: 
g’ (2) =6r° + 3,215 2 — 10,476 2? +9 2? — 0,126 z — 1,101 


6 3,215 — 10,476 9 — 0,126  — 1,101 
— 4,0024 — 5,5787 4,8716 1,8976 
5,148 7,1754 — 6,2660 — 2,4407 g (z,) = 2,3049 z, + 0,7966 
6 8,363 — 7,3030 — 2,8447 2,3049 0,7966 — 1,7854 + 1,6019. 
Als bessere Näherung für die Wurzel erhält man 
9 
— = 0,4256 + 0,6948. 
Karlsruhe. L. Collatz. 174 
Einige Integrale der Besseischen Funk- zliedweise integriert, ergibt sich: 
tionen reellen positiven Argumentes. 
gleichmäßig konvergente Reihe der Funktion): 
und 
z\v 
1 (5) x \2k z\v 
Jy(z) 2k+1 4k 


1) Vgl.E.Brixy: Logarithmus der Besselschen Funk- 
tionen reellen positiven Argumentes. Z. angew. Math. Mech. 
Bd. 19 (1939), S. 373, Formel (5). 


Für v=0 und v=1 erhält man speziell: 
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dr a 
"Ak 1.101140, 
0 
1 1 
== — 1.418865 , 
0 0 


auf 6 Dezimalstellen genau. 


Jv-ı (x) 

2. Integral > 0) d — Durch glied- 
weise Integration der gleichmäßig konvergenten 
Reihe?): 

(2) 


ergibt sich: 


2 (2) 


- 
Iw-1) \2 2k+1\2 


wo die Beziehungen: 


8), 


berücksichtigt wurden. 
Aus der Reihe (3) folgt: 


1 
ı 
Jv (2) 


und für v=P0: 


Mit Hilfe der Relationen: 


dr 


1 1 
u 


ergibt sich weiterhin: 
1 


—1gJ,(1) + 0,1808860. 


Aus der Formel?): 


2 
2 


errechnet man: 


2) Vgl. die Formel (12) des zitierten Aufsatzes. 
#) Vgl. die Formel (29) des zitierten Aufsatzes. 


lg J, (1) = — 0,267621 


und man erhält: 


fie J, (2) d x = 0,913266 — 1 


1 
1 
dz — 0,086734 


beides auf 6 Dezimalstellen genau. 


Die Reihe (4) gibt für v1: 


auf 6 Dezimalstellen genau. 


Jv + 1 (2) 
3. Integral =) dx.— Durch Integration 
der gleichmäßig konvergenten Reihe®): 


erhält man die Ausdrücke: 
Jv 1 (z) d 


Ile) Jv(2) 
2K+2 (5), 
für |z]<2 


Für v=0 und v=1 erhält man die Werte: 
Jı(z) 


Js(z) $ Cı,k 
Hi 
0 


auf 6 Dezimalstellen genau. 


— 0,279773 , 


0,127718 


4. Integral ig dz. - Die Inie- 
gration der gleichmäßig konvergenten Reihe°): 


für |2]<2 und »>0 
führt zu dem Ausdruck: 


4) Vgl. die Formel (19) des zitierten Aufsatzes. 
5) Vgl. die Formel (28) des zitierten Aufsatzes, 


| 
| (4) 
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für |2]<2, v>0 


Daraus ergibt sich: 


fiel 


Il) J,(2) 2 


und für v=0, v=1 speziell: 


(8) 


18 [J, 0,913266 — 1, 


1 


wo die beiden Werte auf 6 Dezimalstellen genau 
errechnet -wurden. 


Die letzten zwei Integrale kann man auch so 
schreiben : 


f* dz— 0,086734 , 


1 
fie 0.042204 . 
0 


Nach Einsetzung von v—2 in die Reihe (8) er- 
hält man: 


1 
8 
36 
ı 
1 
fı — 0,029517 
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BUCHBESPRECHUNGEN 


Prof. Dr. AUGUST THUM VDI, Darmstadt, und 
Dr.-Ing. HANS-RUDOLF JACOBI, Troisdorf, Me- 
chanische Festigkeit von Phenol- 
Formaldehyd-Kunststoffen. Mitteilung 
der Materialprüfungsanstalt an der Technischen 
Hochschule zu Darmstadt (VDI-Forschungsheft 396, 
Beilage zu „Forschung auf dem Gebiete des In- 
genieurwesens“, Ausgabe B, Bd. 10, Mai/Juni 1939). 
39 S. m. 136 Bild. u. 8 Zahlentafeln. Berlin 1939, 
VDI-Verlag G.m.b.H. Preis brosch. 5 M. 


Die Herstellung der von der chemischen Industrie 
erzeugten Kunststoffe, die man seit etwa 80 Jahren 
kennt und die heute unter unzähligen Handels- 
namen in Gebrauch sind, hat in der letzten Zeit eine 
geradezu stürmische Entwicklung genommen. Schon 
vor dem Weltkrieg wurden sie zu allerlei nützlichen 
Gegenständen verwendet und in der Elektrotechnik 
hatten sie sich infolge ihres Isoliervermögens be- 
reits einen Platz gesichert. Nach dem Weltkrieg 
aber und vor allem durch den Vierjahresplan ge- 
lang es dem Chemiker und Ingenieur derartig viele 
Vorzüge dieser Kunststoffe herauszuarbeiten, daß sie 
zu einem Werkstoff von wertvollen mechanischen 
und chemischen Eigenschaften geworden sind. Es 
ist daher nicht verwunderlich, daß man sie nun- 
mehr in allen Zweigen der Technik antrifft, sei es 
in der chemischen oder Maschinenindustrie, in der 
Elektro- oder Bautechnik oder sei es in den 
zahlreichen Anwendungsmöglichkeiten des Haus- 
halts und des täglichen Lebens. Vor 10 Jahren 
war man bestimmt noch nicht überzeugt davon, 
daß mit Lagerschalen und Zahnrädern aus Kunst- 
stoffen namhafte Leistungen übertragen werden 
können. Die Folge dieser Errungenschaften auf 
dem Gebiet der Kunststoffe sind einmal ihre Typi- 
sierung und Normung sowie die Einführung anderer 
z. B. hygienischer Bestimmungen und dann die 
Aufstellung von Prüfverfahren, die die mechanische, 
elektrische, thermische, optische und chemische 
Untersuchung umfassen. 


Eine dieser Untersuchungen, nämlich die mecha- 
nische Festigkeit von Phenol-Formaldehyd-Kunst- 
stoffen haben nun Prof. Thum und Dr.-Ing. Jacobi 
durchgeführt und teilen die Ergebnisse in dem VDI- 
Forschungsheft 396 mit. In den Hauptabschnitten 
dieser Arbeit, die 136 Abbildungen und 8 Zahlen- 
tafeln enthält, wurde zunächst an unbewehrten Stä- 
ben die Festigkeit bei ruhender, bei zügiger und 
bei wechselnder Belastung ermittelt und dann an 
bewehrten Proben, insbesondere die Biegefestigkeit 
im statischen und im Dauerzustand geprüft. Bereits 
aus diesen wenigen Worten, mit denen der Arbeits- 
plan umrissen ist, geht hervor, daß die Verfasser 
ihre Aufgabe in erschöpfender Weise in Angriff ge- 
nommen und auch gelöst haben. Unter dem zahl- 
reichen Schrifttum, das sie mühevoll zusammenge- 
stellt und beigefügt haben, befinden sich allein 
14 eigene Versuchsarbeiten, die Prof. Thum und 
seine Mitarbeiter im Laufe der letzten Jahre durch- 
geführt haben. In Anbetracht der reichen Erfah- 
rungen und tiefen wissenschaftlichen Erkenntnisse 
offenbart sich daher in jedem Abschnitt des vor- 
liegenden Heftes die volle Überlegenheit, mit denen 
seine Verfasser an die Versuche herangetreten sind. 
Die Fachwelt wird das Heft 396 freudig begrüßen, 
und jedem Ingenieur, der gewillt ist, einen neuzeit- 
lichen Bau- und Werkstoff in seinen besonderen 
Festigkeiten kennenzulernen, kann das Studium 
dieser Forschungsarbeit nur empfohlen werden. 

Dresden. Findeisen. 197 


Dr. se. techn. ALFRED MANGER, Der durch- 
laufende Balken auf elastisch dreh- 
barenundelastischsenkbaren Stützen 
einschließlich des Balkens auf stetiger 
elastischer Unterlage. (Mitteilungen aus 
dem Institut für Baustatik a. d. Eidg. Technischen 
Hochschule in Zürich, herausgegeben von den Pro- 
fessoren Dr. M. Ritter und Dr. F. Stüssi, Mit- 
teilung 10.) 170 S. m. 77 Fig., 25 Tab., 8 graph. 
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Beilagen. Zürich und Leipzig 1939, Verlag A.-G. 
Gebr. Leemann & Co. Preis brosch. 7,50 M. 

Die statische Untersuchung eines durchlaufenden 
Trägers mit elastisch drehbaren und elastisch senk- 
baren Stützen ist ein bekanntes Problem der Bau- 
statik, das mehrfach in der Literatur behandelt und 
gelöst worden ist, um die Rechnung dem werk- 
tätigen Ingenieur zu erleichtern und damit die Er- 
gebnisse der Theorie bei der wirtschaftlichen Ge- 
staltung der Tragwerke auszunutzen. Ob jedoch 
dabei viel im Vergleich zur systematischen Lösung 
der Aufgabe nach den allgemeinen Ansätzen zur 
Berechnung hochgradig unbestimmter Tragwerke 
gewonnen wird, ist zweifelhaft. Der Verfasser hat 
sich bemüht, seine Untersuchungen auf alle mit 
dem Problem zusammenhängenden baustatischen 
Fragen auszudehnen, also auch Sonderfälle be- 
handelt und eine Berechnung des Trägers mit 
stetiger elastischer Stützung bei veränderlicher 
Bodenkonstante angegeben. Die theoretischen Be- 
trachtungen werden durch Tabellen, Zahlenbeispiele 
und zahlreiche sorgfältige Zeichnungen ergänzt. 
Die Arbeit wird damit zu einem wertvollen Beitrag 
zur graphischen Untersuchung dieser hochgradig 
statisch unbestimmten Konstruktionen. 


Dresden. K. Beyer. 19 


Dr.-Ing. habil. BERNHARD FRITZ, Dozent a. d. 
Techn. Hochschule Karlsruhe, Bereehnung von 
Kreisgewölben. Gebrauchsfertige Ta- 
bellen zur statischen Berechnung 
von Gewölben mit kreisbogenförmi- 
ger Achse. IV+34 S. m. 12 Textabb. Berlin 
1940, Verlag Wilhelm Ernst & Sohn. Preis 3,20 M. 


Die Arbeit des Verfassers besteht in der. Aus- 
wertung der bekannten Ansätze zur Berechnung 
der drei statisch unbestimmten Größen eines beider- 
seits starr eingespannten symmetrischen Stab- 
bogens mit gleichbleibender Krümmung und gleich- 
bleibendem Querschnitt für ausgezeichnete Stel- 
lungen der Einzellast, ohne jedoch die Symmetrie 
des Tragwerks auszunutzen. Die notwendigen Inte- 
grationen werden dabei numerisch durchgeführt, 
obgleich die strengen Lösungen in der Literatur 
bekannt sind und verwendet werden. Die Tabellen 
werden für 5 verschiedene Pfeilverhältnisse des 
Stabbogens berechnet und enthalten neben den 
Einflußordinaten für die drei statisch unbestimmten 
Größen auch diejenigen für die Biegemomente am 
Bogenkämpfer, Bogenviertel und Bogenscheitel. 
Ihre Anwendung wird mit der statischen Unter- 
suchung von Bogenbrücken gezeigt und an Hand 
mehrerer Beispiele erläutert. Die Arbeit ist daher 
für die Arbeiten des werktätigen Bauingenieurs 
recht brauchbar. 


Dresden. K. Beyer. 19 


Dr. phil. HEINRICH HECHT, Schalt- 
schemata und Differentiaigleichun- 
en elektrischer und mechanischer 
Schwingungsgebilde VI+ 125 S. m. 
51 Abb. Leipzig 1939, Verlag Johann Ambrosius 
Barth. Preis brosch. 9,60 M. 


Wer — vor allem als Leiter einer größeren Zahl 
von Mitarbeitern — sich lange Zeit hindurch mit 
Fragen eines technischen Sondergebietes beschäf- 
tigen muß, wird bestrebt sein, die auf dieses Son- 
dergebiet bezüglichen Sätze und Regeln, die man 
sich zur Lösung der technischen Aufgaben immer 
wieder vor Augen halten muß, zu formulieren; ja 
er wird diesen Sätzen und Regeln unter Um- 
ständen auch besonders angepaßte Fassungen 
geben. Auf diese Weise entstehen die allenthalben 
gebrauchten und nützlichen „Arbeitsblätter“ mit 
ihren Hinweisen, Faustregeln und „Hausgesetzen“. 
Es liegt in der Natur der Entstehung solcher Hilfs- 


mittel, daß sie meist keinen Anspruch auf allge- 
mein gültige Formulierung machen und machen 
können. 

Die vorliegende Schrift ist aus Überlegungen der 
oben genannten Art entstanden. Unter Hinweis auf 
Irrwege, aufgetretene Paradoxa und deren Auf- 
lösung erzählt sie im Plauderton von den physi- 
kalischen Gedankengängen, zu denen die tech- 
nischen Probleme des Baues akustischer Geräte 
Veranlassung gaben. Der Betrachtung unterworfen 
werden verbundene elektrische Schwingungskreise 
und ihre mechanischen Analoga. d. s. Systeme von 
zwei Freiheitsgraden vom Typ der sogenannten 
Ketten. Die Fragen betreffen die Reihen- und 
Parallelschaltung der Elemente und die Art der 
auftretenden Kopplung. Die Behandlung ist voll- 
ständig auf die gerade vorliegenden Verhältnisse 
zugeschnitten. So wird z. B. bei der Frage nach 
der Art der Kopplung abgelehnt, die Entscheidung 
aus den Differentialgleichungen zu treffen; der 
Kopplungsart wird vielmehr physikalische Realität 
zugeschrieben. Dieser Standpunkt läßt sich unter 
bestimmten Voraussetzungen bei der Betrachtung 
von „Ketten“ allenfalls einnehmen, er kann aber 
niemals Anspruch auf alleemeine Gültigkeit erheben. 


Berlin-Charlottenburg. K.Klotter VDI. 167 


RICHARD DOERFLING, Mathematik für 
Ingenieure und Techniker 53 S. m. 
290 Bildern. München und Berlin 1939, Verlag R. 
Oldenbourg. Preis geb. 9,60 M. 


Das Buch bringt dem Ingenieur und Techniker 
die für seinen Beruf unentbehrlichen mathema- 
tischen Grundlagen; es werden die wichtigsten 
Stoffe der Arithmetik und Algebra, Elementar- 
geometrie, Trigonometrie, analytischen Geometrie, 
Differential- und Integralrechnung, Differential- 
gleichungen und Vektoranalysis behandelt. An die 
Vorbildung des Lernenden werden nur geringe An- 
sprüche gestellt; sie beschränken sich etwa auf die 

enntnisse der auf der Mittelstufe einer höheren 
Lehranstalt behandelten Stoffe: diese elementaren 
Dinge werden zudem am Anfange der einzelnen 
Kapitel noch einmal knapp zusammengestellt. 
F'erner hat Verf. sich Leser gedacht, die keine be- 
sondere Veranlagung für Mathematik haben, die 
aber die Mathematik in ihrem Berufe nicht ent- 
behren können, so daß sie sie lediglich ihrer An- 
wendungen wegen studieren wollen. Diese Vor- 
aussetzung macht die in dem Buche eingeschlagene 
Behandlungsweise begreiflich: da es sich in erster 
Linie darum handelte, den Stoff den Lesern ver- 
ständlich und schmackhaft zu machen, war es nicht 
nötig, äußerste Strenge der Beweisführung oder 
lückenlosen systematischen Aufbau anzustreben. 
Wohl aber findet man in dem Buche alles für den 
Techniker Wichtige einfach und klar zur Dar- 
stellung gebracht, insbesondere trifft man auch 
Stoffe an, die sonst in solchen elementaren Ein- 
führungen nicht immer behandelt werden, deren 
Kenntnis aber für den Techniker unerläßlich ist, 
so die Lehre von den Hyperbelfunktionen, von den 
Fourierschen Reihen, die Theorie der wichtigsten 
gewöhnlichen und partiellen Differentialgleichungen, 
den Stokesschen Satz usw. Zur Festigung und Er- 
läuterung des mathematischen Stoffes werden zahl- 
reiche Anwendungen gebracht, u. a. findet man 
eine kurze Einleitung in die Mechanik, die dann in 
der Lehre von den Differentialgleichungen An- 
wendungsbeispiele liefert, ferner Anwendungen der 
Vektoranalysis auf die Lehre von der Gravitation 
und die Theorie des elektromagnetischen Feldes 
u. a.m. Das Buch dürfte eine recht geeignete Ein- 
führung in die Ingenieurmathematik darstellen, für 
das Selbststudium wie zum Gebrauch neben den 
Vorlesungen geeignet. 


Berlin. E. Mosch. 168 
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Dr. K. SCHWIDEFSKY, wissenschaftl. Mitarbei- 
ter d. Opt. Werke Zeiß, Jena, Einführungin 
die Luft- und Erdbildmessung. 2. erw. 
u. verb. Aufl. m. 73 Abb., drei schwarzen und zwei 
farbigen Tafeln im Text, einer schwarzen Tafel, 
einer farbigen Brille und zwei Stereobildern im An- 
hang. IV + 137 8. Leipzig und Berlin 1939, Verlag 
B. G. Teubner. Preis geb. 8 M. 


Die Bildmessung ist bekanntlich ein topographi 
sches Verfahren, dessen hervorragende Eignung für 
zivile, militärische und koloniale Vermessungen 
allgemein erkannt ist. Darüber hinaus aber löst die 
Bildwessung in der messenden Wissenschaft und 
Technik wichtige Aufgaben, so z. B. im technischen 
Versuchswesen, in der Ballistik, Physik, Röntgen- 
medizin und Tierzucht. 

Das vorliegende Buch, dessen erste Auflage an 
dieser Stelle (Bd. 18 [1938], S. 259) besprochen 
wurde, ist in seiner rasch notwendig gewordenen 
zweiten Auflage durch wesentliche Ergänzungen 
auf den neuesten Stand der Technik gebracht wor- 
den. Mustergültig in seiner knappen, wissenschaft- 
lich einwandfreien und dabei verständlichen Dar- 
stellung, mustergültig auch in seinen instruktiven 
Abbildungen, ist das Buch die beste Einführung in 
dieses in seinen mathematischen, optischen, mecha- 
nischen und physiologischen Grundlagen gleich 
interessante Gebiet. 


Dresden. Hugershoff. 182 


Dr. RICHARD FINSTERWALDER, a. o. Prof. a. 
d. Techn. Hochschule Hannover, Photogram- 
metrie. 237 S. m. 103 Abb. u. 17 Tabellen. 
Verlag Walter de Gruyter &Co. Preis 
geb. 14 M. 


Das vorliegende Lehrbuch nimmt, entsprechend 
der Absicht des Verfassers, eine Mittelstellung ein 
zwischen den bekannten umfangreichen Hand- 
büchern von Bäschlin, Gast, v. Gruber und Hugers- 
hoff einerseits und der ausgezeichneten Einführung 
von Schwidefsky andrerseitt. — Das Erscheinen 
dieses neuen Lehrbuches ist sehr zu begrüßen, weil 
die oben erwähnten umfassenden Handbücher in- 
folge der überaus raschen Weiterentwicklung der 
Bildmessung hinsichtlich der Arbeitsverfahren schon 
nicht mehr ganz dem heutigen Stand der Technik 
entsprechen. Das wohlgelungene Buch behandelt 
deshalb auch besonders eingehend die Haupt- 
aufgabe der Photogrammetrie, 
nämlich die Doppelpunkteinschaltung im Raum, 
ferner die Aero- und Radialtriangulation, die für 
militärische und koloniale Aufgaben wichtigen 
Näherungs- und Schnellverfahren und die graphische 
und optisch-mechanische Entzerrung. Die mathe- 
matischen Grundlagen für diese Verfahren und die 
in ihnen wirkenden Fehlerquellen sind unter Ver- 
wendung sehr geschickt entworfener Abbildungen 
leicht verständlich dargestellt. 
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Ferner sind bei der Schriftleit 
Bücher eingegangen (ausführliche 
bleibt vorbehalten): 


Dr. phil. WILHELM SPA Physik und 
Technik der Härte und Weiche. VII+ 
250 S. m. 214 Textabb. Berlin 1940, Verlag Julius 
Springer. Preis geb. 19,50 M. 


Prof. Dr. A. WITTING, Integralrechnun 
(Sammlung Göschen, Bd. 88). 2. verb. Aufl. 175 S. 
m. 62 Fig. u.:190 Beispielen. Berlin 1940, Verlag 
Walter de Gruyter & Co. Preis geb. 1,62 M. 


folgende 
prechung 


Dr. JOHANNES TROPFKE, Geschichte der 
Elementar-Mathematik insystemati- 
scher Darstellung mit besonderer 
Berücksichtigung der Fachwörter 
IV. Bd. Ebene Geometrie. 3. verb. u. vermehrte 
Aufl. besorgt v. Dr. Kurt Vogel. 316 S. Berlin 
1940, Verlag Walter de Gruyter & Co. Preis 
geb. 13 M. 


Dr.-Ing. E. h. WOLMAR FELLENIUS, Prof. a. 
d. Kgl. Techn. Hochsch. zu Stockholm, Erdsta- 
tischeBerechnungenmitReibung und 
Kohäsion (Adhäsion) und unter An- 
nahme kreiszylindrischer Gleit- 
flächen. 2. ergänzte Aufl. IV +48 S. m. 38 Text- 
abb. Berlin 1940, Verlag Wilhelm Ernst & Sohn. 
Preis brosch. 4,80 M. 


Prof. Dr. HERMANN SCHUBERT, Mathema- 
tische Mußestunden, eine Sammlung 
von Geduldspielen, Kunststücken 
und Unterhaltungsaufgaben mathe- 
matischer Natur. 6. Aufl, neubearb. von 
Prof. Dr. F. Fitting. 260 S. Berlin 1940, Verlag 
Walter de Gruyter & Co. Preis geb. 4,80 M. 


Dr. phil. GERHARD HESSENBERG, Ebene 
und sphärische Trigonometrie. 4. Aufl. 
(Sammlung Göschen, Bd. 99.) 171 S. m. 59 Fig. 
Berlin 1940, Verlag Walter de Gruyter & Co. 
Preis geb. 1,62 M. 


EUGEN DIESEL, Das Phänomen der 
Technik, Zeugnisse, Deutung und 
Wirklichkeit. 2. Aufl. 259 S. m. 35 Bildtafeln. 
Berlin 1940, VDI-Verlag, G.m.b.H., und Leipzig, 
Verlag Philipp Reclam jun. Preis geb. 7,50 M. 


Dr. KARL KOMMERELL, em. o. Prof. a. d. Uni- 
versität Tübingen, Vorlesungen über ana- 
Iytische Geometriedes Raumes. VII+ 
388 S. m. 77 ji © im Text u. auf 2 Tafeln. Leipzig 
1940, K. F. Koehler Verlag. Preis geb. 20 M. 


HEINRICH DÖRRIE, Grundrißder Physik 
mit besonderer Berücksichtigung der 
Anwendungen. XII + 836 S. m. 456 Fig. Bres- 
lau 1940, Verlag Ferdinand Hirt. Preis geb. 20 M. 


Dr.-Ing. HEINRICH PRESS, Der Boden als 
Baugrund mit Ergebnissen eigener 
Versuche. 2. neubearb. Aufl. (Mitteilungen 
aus dem Gebiete des Wasserbaues und der Bau- 
grundforschung, Heft 9. V-+49 S. m. 64 Text- 
abb. Berlin 1940, Verlag Wilhelm Ernst & Sohn. 
Preis brosch. 5,80 M. 


Dr. FR. REGER, Tachymetertafeln für 
neue (zentesimale) Teilung. + 
490 S. Stuttgart 1940, I. B. Metzlersche Verlags- 
buchhandlung. Preis geb. 19,40 M. 


Prof. WILHELM SPANNHAKE und Dipl.-Ing. 
HERMANN DOSENBACH, Geschwindigkeit 
und DruekinSchaufelkanälen hydro- 
Getriebe. 1. Teil. (Deut- 
sche Kraftfahrtforschung, Heft 38.) 16 S. m. 26 Bil- 
dern. Berlin 1940, VDI-Verlag G.m.b.H. Preis 
brosch. 1,40 M. 


Prof. Dr. ERNST RAUB, Dipl.-Ing. GERHARD 
BUSS u. Dipl.-Ing. MAX ENGEL, Die Verchro- 
mung von Fahrzeugzubehörteilen, 
Austausch der Vernicklung gegen 
die Verkupferung vor der Verchro- 
mung. (Deutsche Kraftfahrtforschung, Heft 40.) 
29 S., 31 S. m. 30 Bildern. Berlin 1940, VDI-Ver- 
lag G.m.b.H. Preis brosch. 3 M. 
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